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Abstract 

We study Yukawa unification, including the effects of a pliysical neutrino mass consis- 
tent with the Superkamiokande observations, in a string/D-brane inspired S'f/(4) 
SU{2)i ® SU{2)fi model which allows the most general non-universal scalar and 
gaugino masses, including the usual D-term contributions which arise in 5*0(10). We 
investigate how the tight constraints from rare decays such as 6 — > S7 and t —>■ fij 
can provide information about the family dependent supersymmetry breaking soft 
Lagrangian, for example the trilinears associated with the second and third family. 
Many of our results also apply to SO (10) to which the model approximately reduces 
in a limiting case. In both models we find that Yukawa unification is perfectly viable 
providing the non-universal soft masses have particular patterns. In this sense Yukawa 
unification acts as a window into the soft supersymmetry breaking Lagrangian. 
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I. INTRODUCTION 



One of the earliest successes of the Minimal Supersymmetric Standard Model 
(MSSM) was that it could allow for the unification of the gauge couplings at 
high energy {Q = Mx) 0, thereby opening the door for Supersymmetric Grand 
Unified Theories (SUSY GUTs) § |, |, §. SUSY GUTs typically involve some third 
family Yukawa unification 0, ^ |10|. For example minimal SU{5) predicts 
bottom-tau Yukawa unification at Mx (A;, = A,-) and allows for the prediction of the 
top quark mass (m^) in terms of the ratio nih/mr. In models like 5*0(10), and in the 
SU{A)®SU{2)l®SU{2)r (422) Pati-Salam model 0, complete third family Yukawa 
unification holds (Ai = Af, = Ar). Thus the ratio of the vacuum expectation values 
(VEVs) of the up/down Higgs bosons doublets - tan/? = ^2/^1 - can be predicted. 
Since the top and bottom Yukawa couplings start out equal at the GUT scale, and 
their renormalization group (RG) evolution is very similar (they differ only by the 
f/(l) couplings at one loop) their low energy values are approximately equal so we 
expect tan/3 ~ mt/rrib ~ 40. 

In this paper we shall include the neutrino Yukawa coupling in the analysis giving 



quadruple Yukawa unification (A^ = A;, = A^ = X^) as was first done in Ref. |10|. In 
the light of the Superkamiokande results [|13| (assuming hierarchical neutrino masses 
and the see-saw mechanism) we regard the third family neutrino mass {m,^) as an 
input which allows for an additional prediction - the mass scale at which the heavy 
right-handed neutrinos decouple - M,y. However such predictions depend sensitively 
on the parameters of the soft SUSY breaking Lagrangian, as we now discuss. 

Yukawa unification (assuming the MSSM and SUSY GUTs) at first was thought 
to lead to an acceptable top mass prediction, given the experimentally permitted 
range of top and bottom quark masses and a^, with a unified Yukawa coupling at 
Mx of order one (Ax ~ 1) @, 0' 13' EH HI' Ol- However, when low energy SUSY 
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corrections to the running bottom quark mass due to the decouphng of SUSY particles 
were included []T8|, ^ , the initial good agreement was spoiled. The SUSY corrections 
to the bottom quark mass are : 

6mh /i tan (3 



Arrih 



o \2 

r/ 2 2 2 \ , A T/^ 2 2 2 



(1) 



where the function I{x,y,z) is given in [18|. The corrections can be positive or 



negative depending on the sign of fi, and being proportional to tan/5, are rather 
large. 

With universal soft masses and negative fi, the SUSY corrections to the bottom 
mass lead to large negative Ami, = 5mi,/mi, ~ —20% corrections which implies a 
larger h quark mass before the SUSY corrections, which effectively lowers the previ- 
ously successful top mass prediction to ~ 150 GeV which is too small |TB|. With 
universal soft masses and positive ^ the top mass is predicted to be too large, outside 
its perturbative upper limit. It has already been pointed out in the literature that 
5mfe can be made small by either assuming explicit non-universal scalar soft masses 



21| , or by introducing approximate Peccei-Quinn and R symmetries ^3 
However in the framework of 50 (10) the degree of non-universality one could assume 
until recently appeared rather limited. For example, if both Higgs doublets arise from 
a single 10 dimensional representation then they will necessarily have a common soft 
scalar mass at the GUT scale, and the same applies for the light Higgs bidoublet of 
the Pati-Salam model. Such universality is also a problem for electroweak symmetry 
breaking where one requires a large hierarchy of vacuum expectation values (VEVs) 
starting from very symmetrical initial conditions where the two Higgs doublet soft 
masses are equal, and where the approximately equal top and bottom Yukawa cou- 
plings tend to drive both Higgs masses negative equally, making large tan/3 rather 
difficult to achieve. 

A large step forward for both these problems has been to realise the importance 



of D-teim contributions to scalar masses EM, which naturally split the two Higgs 



doublet masses. If the up-type Higgs doublet mass at the GUT scale is smaller than 
the down-type Higgs doublet mass, then this makes electroweak symmetry breaking 
with large tan (3 much more natural. Assuming negative /i, the D-terms also allow 
a choice of non-universal scalar masses which reduces the correction to the 6-quark 
mass and hence allow a larger top quark mass. Although these problems appear to 
be resolved with D-terms, one still faces difficulties with rare decays such as b ^ 57, 
which is also enhanced for large tan/3, and negative ^ [^. This problem can be 
avoided either by increasing the masses of all the superpartners, or by considering 
positive fi which will tend to cancel the SUSY contributions, however both these 
procedures lead to fine-tuning. Another possibility which we pursue in this paper 
is to consider the effect of non-universality in the family space, which may lead to 
additional contributions to 6 ^ 57 which can cancel those coming from the family 
universal sources. We explicitly check that the contributions that we introduce do 
not introduce problems elsewhere such as with r — »• fi^. 

The main purpose of this paper is to provide a detailed analysis of Yukawa uni- 
fication, post-Superkamiokande, allowing the most general non-universal soft SUSY 
breaking masses possible. Apart from non-universality in the soft scalar masses aris- 
ing from D-terms, we shall also consider more general types of non-universality which 
may arise in models such as the string/D-brane inspired 5'f/(4) (g) SU(2)i ® SU(2)r 



model [^],|^, which reduces to the 5*0(10) model in a limiting case. This allows 
more general soft scalar mass non-universality including similar D-term contributions 
as in S'O(IO) (since the relevant broken f/(l) factors all arise from the breaking of the 
SU {A)®SU{2)l®SU{2)r subgroup), and also permits violations of gaugino mass uni- 
versality which have not so far been studied. The underlying theme of our approach 
is that the sensitivity of Yukawa unification to soft SUSY breaking parameters is to 
be welcomed, since it provides a window into the soft SUSY breaking Lagrangian, 
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both from the point of view of family universal and family non-universal soft SUSY 
breaking parameters. 

The layout of the remainder of the paper is as follows. In the section II we present 
our calculational approach, and in section III we review the SUSY corrections to 
the bottom and tau masses. Section IV deals with the effect of neutrino Yukawa 
couplings on Yukawa unification and in section V we present results for gauge and 
Yukawa unification, assuming minimal supergravity (mSUGRA) with universal soft 
SUSY breaking parameters and including the neutrino threshold in addition to the 
usual low energy thresholds. In section VI we turn to Yukawa unification in the 
SU{4:) (g) SU{2)l ® SU{2)r model, including the effect of physics above the GUT 
scale, and non-universal soft scalar and gaugino masses permitted in this model. We 
also consider the effect of D-terms in the 5'O(10) hmit of the model, and then study 
b ^ sj and t —>■ /ij including non-universality in the family- dependent trilinear 
parameters. Section VII concludes the paper. 
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II. FRAMEWORK 



In this section we will summarize how we implemented the RGEs and matching 
boundary conditions of the model. In the region Q < Mz the effective theory is 
5f/(3)c ® t/(l)em, thus 3-loop QCD [0 plus 2-loop QED RGEs apply. Between 
Q = Mz and Q = Ms, where Ms is the scale that parameterises the energy at 
which the theory effectively becomes supersymmetric (see Appendix B for details), 
we considered 2-loop Standard Model (SM) RGEs in the gauge and Yukawa couplings 
H, m. In the region Ms < Q < Mx we evolved 2-loop (1-loop) gauge/Yukawa 



(all other parameters) MSSM RGEs [|T^, 31|, properly adapted and extended to take 
into account the presence (and decoupling) of right handed neutrinos z/^. 

A. Input 

The low energy input was : the running electromagnetic coupling a~^{Mz) = 
127.8, the pole tau lepton mass M^- = 1.784 GeV and several values for the pole 
bottom quark mass M^ = 4.7, 4.8, .., 5.1 GeV. We have converted the pole bottom 
mass into "running" bottom mass in the MS scheme using two loop QCD perturbation 
theory : 



m,(g = M,) = ^ (2) 

^ , 4 a,(M,) , fasjM,)-^ 

i + h Kh 



3 TT \ 71 

where Kf, = 12.4. The corresponding running masses, for = 0.120, are mh{Mii) = 
4.06, 4.15, 4.24, 4.33, 4.42 GeV. The experimental range for estimated from bot- 
tomonium and B masses is 4.0-4.4 GeV The strong coupling as{Mz) was taken 
to be in the range 0.110-0.130. 

The low energy input was complemented by the Super-Kamiokande atmospheric 
neutrino data that suggests ml — ml = 10^^ to 10^^ eV^ |13]. Assuming that the 



1^3 1^2 

neutrino masses are hierarchical we required that m^.^ ~ 0.05 eV. 



The universal high energy inputs assumed initially are motivated by mSUGRA: 
spontaneously broken supergravity in which the local N=l super symmetry breaking 
occurs in a "hidden" sector and is only transmitted to the "visible" sector through 
weak gravitational - flavour blind - interactions. Thus we assume to begin with 
universal soft SUSY breaking masses (USBM) given by a common gaugino mass 
Ml/2 and soft scalar masses ttiq a.t Q = Mx- The trilinear A-terms were set to zero. 
In section VI we relax the universality assumption. 

In the Yukawa sector we assumed that the third family Yukawa couplings have 
their origin in a unified renormalizable interaction which fixes their values to be Ax- 

B. Running and matching boundary conditions 

The process by which the output was generated relied in the initial estimation 
and successive iterative refinement of a priori unknown parameters such as Mt and 



tan/?. The procedure is described in detail in section VII of Ref. |jT6|. 

Starting at Q = Mz, as{Mz) and ae{Mz) were firstly run down to Q = 1 GeV 
(where the up, down and strange quark running masses were fixed |]29|) and secondly 
up to Q = Mz using SU{3)c ® U{l)em RGEs. In the "running up" process the heavy 
(charm, tau and bottom) fermion pole masses were converted to running ones using 
the expressions in Ref. [^. At Q = Mz, o.^ was corrected for the decoupling of heavy 
gauge bosons, Higgs and Nambu-Goldstone bosons, and ghosts ||33|. Afterwards the 
value of si |^ was used to obtain g' and g. The bottom and tau Yukawa couplings, in 
the Standard Model, were evaluated using \^^{Mz) = mi,^r{Mz)/v, where the VEV 
is f = 174 GeV. Next, the gauge and the Yukawa couplings were run from Q = Mz to 
the estimated pole top mass Q = Mt using 2-loop SM RGEs, at which point Mt was 
converted to MS running mass mt{Mt), joining the list of parameters to be integrated 
to Q = Ms- All along, threshold corrections in the gauge couplings were included by 
changing the 1-loop /^-functions (using the "step" approximation.) At Q = Ms the 



MS gauge couplings at were converted to the DR ones 0, ||, Q and corrected for the 



cumulative effect of decoupling of all the SUSY particles [0. The top, bottom and 
tau Yukawa couplings were then converted from the SM to the MSSM normalization 
and corrected for the SUSY corrections : 



^MSSM ^ Af^/sin/? (3) 
^MSSM ^ AfVcos/3-5m,/^i (4) 
^MSSA4 ^ AfVcos/3-5m,M (5) 

Afterwards we run up all the above couplings, together with an estimate for 
X^^'^^^^Ms), in two stages, firstly from Q = Ms to Q = M^, and afterwards from 
Q = Ml, to Q = Mx, properly excluding (including) z/^ in the former (latter) stage. 
In our model Mx was fixed to be the scale at which only the U{1) and SU{2) gauge 
couplings unify : ai{Mx) = a2{Mx) 7^ asi^Mx)- 

At Mx gauge and Yukawa unification were tested leading to eventual wiser choices 
to the next estimates for Mt, tan/?, Mx and M,^. 

The iteration cycle was completed by setting the USBM to their unification values 
after which all the couplings and masses of the model were run down from Q = Mx 
to Q = Ms using the inverse "running up" procedure described above. Finally at 



Q = Ms the 1-loop effective Higgs potential was minimized |l35| , p6| and the SUSY 
Higgs mixing parameter /i^ and the corresponding soft term mg were determined using 



m, 37, 38 



9 m? — motan^/3 1 ,^9 



-3 = (/^2+/^l)^^^^^^^ (7) 

where /i| ^ = m| ^ + fi'^ and the up/down soft Higgs boson masses are 1 = mfl'^^ + 
^The S2,i parameterise the 1-loop corrections to the tree level Higgs potential [p9[. 



We note that the number of independent parameters the model can predict is 
four - as many as the constraints imposed (one gauge and three Yukawa unification 
conditions.) We took them to be M^, M^, tan/3 and My. The latter was fixed by 
requiring that rriy.^ = 0.05 eV. 

III. SUPERSYMMETRIC BOTTOM 
AND TAU MASS CORRECTIONS 



In the Standard Model it is frequent to work in off-shell mass schemes in which 
the running masses mf{Q) differ from their physical masses Mj (defined as the real 
part of the complex pole position of its propagator) by some finite correction For 
the quarks the most important corrections (arising from gluon loops) are well known 
and particularly affect bottom quark mass (see Eq. @).) In supersymmetric models 
additional corrections are also present. For large values of tan/3 some of these SUSY 



corrections can indeed affect the running bottom mass mj, by 20 % [jl8|, |T9[ , thus their 
consideration is crucial for the prediction of M^. In this section we review the origin 
of the SUSY corrections to the bottom quark (5mfc) and tau lepton mass ((5m^.) 
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FIG. 1 Gluino and higgsino loop diagrams corresponding to the bottom mass cor- 
rections S^mi, and 5^ mi, of Eq. (||) and Eq. (||). 

^L^Tl ^R 

^— i > 

'Cl B 'Cr 

FIG. 2 Bino diagram corresponding to the tan mass correction S^tUt- of Eq. (p^. 

In the MSSM, at tree level, the bottom quark acquires a mass exclusively through 
the non- vanishing VEV of = Vi. However, at 1-loop level, the bottom quark also 
receives a "small" mass from (h^) = V2. The dominant processes responsible for the 
additional corrections to the bottom (tau) mass are characterized by gluino/sbottom 
and higgsino/stop (bino/stau) diagrams which we illustrate in Fig. 1 (and Fig. 2.) 
The corresponding expressions for S^rrii,, 6^mb and 6^171^ are given by [^, [T^, 5T| : 



/X'^rrih = = rria u tan p limp,, rrir , mr 



A^mb = ^^ = ^ A/i tan/3/(/i2,m? ,m|) (9) 



A^mr = = — mo /X tan /3 /(m|,m? ,m? ) (10) 

rriT- An a \ 2 

where rrig {rrij^) is the gluino (bino) mass, at = A^/47r, a' = g'"^ /At:, At the top soft 

trilinear term and ^ - the sbottom, stop and stau masses. The function / is positive, 

symmetric, smallest for degenerate masses and approximately scales with the inverse 

of its biggest argument. We also find convenient to define the total absolute and 
^See footnote in appendix C. 



relative bottom corrections : 

Smb = S'^rub + 5"mb Anib = A^nib + A^mfe. (11) 

The bottom mass before the SUSY corrections are included - rrib = Xb^i - is related 
to the bottom mass after the SUSY corrections are included - mf^^^ - through : 

^susY =rnb + 6mb = mb{l + Amb) (12) 

The pole mass after the SUSY corrections are included is given by Eq. (|) using 
^susY Eq. (|1|). 
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IV. DECOUPLING OF THE HEAVY 
RIGHT HANDED NEUTRINO 



In this section we briefly discuss the decouphng of the heavy SM singlet right- 
handed neutrino with mass M^. We assume that the neutrino Yukawa matrix is 
dominated by a single entry in the 33 position, although in realistic models of neu- 
trino mixing one would expect that there are at two entries which may have similar 
magnitude. We have checked that our results are insensitive to the presence of a large 
off-diagonal entry in the neutrino Yukawa matrix, assuming that the charged lepton 
Yukawa matrix does not have any large off-diagonal entries. Q The decoupling is 
achieved via see-saw mechanism which generates a small mass for the left-handed 
neutrino through the presence of a right-left Dirac Yukawa coupling. The part of the 
superpotential of interest is : 

W = z/" A, z/ /i° + i M, z/^ z/^ (13) 

Thus the light neutrino tau acquires a mass = A^f|/(4M,^). In our model we used 
= 0.05 eV suggested by Ref. O] to fix ~ 10^^ GeV (see results in Table IV 



in section V.C.) 

V. mSUGRA RESULTS 



We now proceed to discuss the results generated by the model described in the 
previous sections. Although many of the results presented here appear elsewhere, we 
find it useful to compile and review them here for the purposes of comparison to the 
new situations we discuss later such as the effect of neutrinos and non-universal soft 
masses. These have been organized in three categories which are suitable to expose 
their variation with as{Mz) in the range 0.110-0.130, selected pole bottom masses 

■^We focus on the third family only, thus we simplify our notation by replacing v-r ^ v and 
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Mb = 4.7, 4.8, .., 5.1 GeV and universal gaugino/soft scalar masses M1/2, mo < 1 TeV. 
We plotted graphs showing the dependence of the results with as for various Mb and 
fixed Ml/2 = 400, uiq = 200 GeV; graphs scanning the Mi/2-mQ parameter space 
for illustrative fixed and M^, and a set of numerical tables corresponding to the 
results obtained from nine models for which the input is listed in Table I. 



TABLE I. 



Case 




Mb 


TTlbiMb) 


Ml/2 


mo 


A 


0.1150 


4.70 


4.12 


400 


200 


B 


0.1150 


4.70 


4.12 


800 


400 


C 


0.1150 


5.10 


4.49 


400 


200 


D 


0.1150 


5.10 


4.49 


800 


400 


E 


0.1250 


4.70 


3.99 


400 


200 


F 


0.1250 


4.70 


3.99 


800 


400 


G 


0.1250 


5.10 


4.36 


400 


200 


H 


0.1250 


5.10 


4.36 


800 


400 


I 


0.1250 


4.80 


4.08 


600 


400 



TABLE I. Input values of as{Mz), the pole bottom mass Mb, the running bottom 
mass mb{Mb) obtained from Mb using two loop QCD corrections only (see Eq. (|^)), 
and the universal gaugino and soft scalar masses M1/2, mo (given in GeV units) at 
the unification scale Q = Mx for a list of models we will refer in the main text as 
Case A,B,..,I. 



A. Gauge unification 



In Fig. 3 we plot the unified U{1) and SU (2) gauge couplings {Mx)=a2^ {Mx) = 
OL^ together with oi^^iMx) against as{Mz) for M1/2 = 400, mo = 200 GeV and sev- 
eral values of Mb- The mismatch Aa^^ = a^^(Mx)~ax^ in gauge unification is small 
and decreasing for increasing a^- In fact, for large values of ~ 0.128 complete gauge 
unification occurs. The sensitiveness of Aa^^^ to Mb is quite small. Complete gauge 
unification can also be present for lower values of ~ 0.125 as is shown by case H 
in Table II where we find ~ 03 ~ 25.18 (M1/2 = 800, mo = 400 GeV.) 
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0.128 



.13 



FIG. 3. Dependence of unified gauge couplings a^^{Mx) = a2^{Mx) = ol^ and 
a'^^{Mx) on as{Mz) for several bottom masses M(, and fixed gaugino and soft 
scalar masses Myii(Mx), mo{Mx)- 



TABLE II. 



Case 






Mx/lO^s 


A 


24.70 


25.46 


1.64 


B 


25.21 


25.83 


1.37 


C 


24.70 


25.38 


1.52 


D 


25.18 


25.75 


1.26 


E 


24.68 


24.87 


1.83 


F 


25.20 


25.25 


1.54 


G 


24.68 


24.80 


1.70 


H 


25.18 


25.17 


1.42 


I 


24.98 


25.09 


1.66 



TABLE II. Predicted values for the unified gauge coupling a^^ = aj^^{Mx) = 
a2^{Mx), for the strong coupling ag ^(Mx) and for the unification scale Mx (given 
in GeV.) 



Analysing Table II we find that Sax^ — ^(^x^/^x^ < 3%. Gauge unification also 
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depends on the SUSY spectrum. Generally, increasing M1/2 and/or mo decreases Mx 
and 5a^. In short partial gauge unification is achieved for ~ 1 — 2 GeV 

with ~ 25.0 ±0.5 while complete unification = a^^{Mx) is favoured by large 
a,{Mz), Ml/2 and M^. 

B. Top mass prediction 
with SUSY bottom and tau mass corrections 

We now turn to the predictions for the top mass. Experimentally the top mass 
has been measured to be |^ : 

Mt ~ 174.3 ± 3.2(stat.) ± 4.0(syst.) GeV (14) 
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'54/ 


y y y y i 


Ml/2 = 400 GeV - 
mo = 200 GeV 

1 



0.11 0.115 0.12 0.125 0.13 

a.(Mz) 

FIG. 4. Predicted values for the pole top mass Mt against as{Mz) for several values 
of the pole bottom mass in the range Mf, = 4.7— 5.1 GeV. The dashed line refers to 
the top mass prediction when no SUSY correction to the bottom mass is considered 
((5mb = 0). The prediction for the top mass when negative corrections to the bottom 
mass are included (/i < 0) is indicated by the solid line {Smi, ^ 0.) 



In Fig. 4 we compare the values for the pole top mass Mt prediction obtained 
when no SUSY correction to the bottom mass is considered {5mi, = 0), indicated 
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with a dashed hne, and when they are included {Smi, ^ and negative) plotted with 
a solid line. We see that the SUSY bottom corrections (for our choice of negative /i) 
decrease the top mass considerably. 

In fact, the magnitude of the SUSY correction is so large that it excludes the 
possibility of an eventual positive sign for fi. The reason is because if we have a pos- 
itive /i the SUSY bottom corrections are positive. Thus, in order to keep the bottom 
mass, after the SUSY corrections are included mf^*^^, in the allowed experimental 
range, the bottom mass before the SUSY corrections are included must decrease. In 
this case, third family unification leads to an increase in the top Yukawa coupling. 
However, it turns out that the required increase in the top Yukawa coupling is so 
large that it drives it, at high energy, well beyond perturbation theory. Numerically, 
we find that the RGEs fail to converge. ^ 

At this point it is also interesting to comment on the effect the bino correction 
to the tau mass has on the prediction for the top mass. Again, a negative /i means 
that S^rrir is negative (see Eq. (p!0D). Thus, in order to keep the tau mass after the 
SUSY correction is included unchanged, the tau mass before the SUSY correction is 
included must increase, implying that, third family unification predicts a larger top 
mass. Numerically we find that the top mass increases by 3-4 GeV. 

Returning to Fig. 4 we observe that when the bottom corrections are included the 
top mass prediction is only acceptable for large values of as{Mz)- For this reason, 
in the analysis that will follow, we will study the implications of including the SUSY 
bottom corrections in our model by taking as{Mz) to be in the range 0.120-0.130. 



^One can turn the argument around by saying that, if the top quark mass is fixed to be around 
175 GeV then third family Yukawa unification at Mx requires, for positive /i, a very large bottom 
mass prediction after the SUSY corrections are included. 
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FIG. 5. Pole top mass prediction Mt against as{Mz) for several values of the pole 
bottom mass M(, and fixed Mi/2) "^o- 




FIG. 6. Prediction for tan/3 against as{Mz) for several values of the pole bottom 
mass Mb and fixed M1/2, mo- 
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In Fig. 5 and 6 the pole top mass Mt and tan (3 are plotted against as for several 
values of the bottom mass M^. We can see that both are significantly sensitive to 
as and to M^. For = 4.7 (the first upper line) the top mass takes values within 
165.6 < Mt < 178.4 GeV and 45.73 < tan/3 < 49.41 for 0.122 < a, < 0.130. 
Additionally for a fixed value of as = 0.125 the top mass increases from 157.7 to 
170.6 GeV and tan/3 from 41.99 to 47.15 when Mf, decreases from 5.1 to 4.7 GeV. In 
Table III we list the numerical predictions for the top mass, tan /? and the value of 
the unified Yukawa coupling Ax- 



TABLE III. 



Case 


Mt 


tan P 


Ax 


A 


152.7 


42.16 


0.426 


B 


158.1 


44.16 


0.469 


C 


135.3 


35.94 


0.317 


D 


139.7 


37.51 


0.340 


E 


170.6 


47.15 


0.569 


F 


176.2 


49.40 


0.643 


G 


157.7 


41.99 


0.428 


H 


163.3 


44.06 


0.473 


I 


170.7 


47.34 


0.561 



TABLE III. Predicted values for the top mass (in GeV), tan/3 and for the value of 
the third family unified Yukawa coupling Ax- 

In Fig. 7 the strong correlation between Mt and Ax is exposed. It is clear that 
since Ax is not large the top quark mass prediction is far from its infrared (IR) fixed 



point - 200 GeV [44, 45 



In Fig. 8 we show that Mt increases with increasing SUSY particle masses. The 
reason for such dependence can be traced to the SU (3) gauge group factors in the SM 
and MSSM RGEs : dAf*'''^''^^''Vdt ~ _bf^^^MSSMg2_ Numerically we have bf'' = 
8 > bf^^^'^ = 16/3 thus increasing Ms allows for a wider Mz < Q < Ms range of 
integration for the SM which favours an enhancement for the top mass. 
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0.5 0.55 0.6 0.65 

Ax 



FIG. 7. Correlation between the top pole mass prediction Mf and the value of the 

unified gauge coupling at Mx, Ax- Each line corresponds to a fixed choice of Mb 
and along it as(Mz) decreases from a maximum value of 0.130 (indicated with a 
crossed circle) to lower values at 0.001 intervals marked with crosses. 




600 605 610 615 620 625 

Ms/GeV 



FIG. 8. Dependence of the top mass prediction Mt on the effective supersymmetry 
scale Ms- Each line corresponds to a fixed choice of M), and along it as{Mz) 
decreases from a maximum value of 0.130 (indicated with a crossed circle) to lower 
values at 0.001 intervals marked with crosses. 
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FIG. 9. Correlation between the top mass prediction Mt and the prediction for 

tan /?. Each hnc corresponds to a fixed choice of toq and along it M1/2 decreases 
from a maximum value of 1 TeV (indicated with a crossed circle) to lower values at 
50 GeV intervals marked with crosses. 

Perhaps more interesting than Fig. 8 where variations in Ms are indirectly induced 
by as is Fig. 9 where as — 0.125 is keep fixed but Mi/ 2 and mo are allowed to vary. 
The results are displayed in the M^-tan P plane. 

A quick estimate of the effect of including a consistent supersymmetric scale Ms 
on the top mass prediction can easily be computed. Taking, for example, mo = 400 
GcV, varying M1/2 from 600 to 700 GeV will increase Mt from 170.7 to 172.0 GeV 
(see Fig. 9) and Ms increases from 915 to 1047 GeV. Thus if we had considered a 
"rigid" Ms ^ Mz ^ 100 GeV we would find a value for Mt decreased by AM* ~ -8 
GeV. 
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C. Decoupling of the heavy tau neutrino 



In Fig. 10 we show an example of third family Yukawa unification generated for 
case I. We note that the inclusion of the SUSY correction 5mi, leads to low values 
for the unified Yukawa coupling Ax ~ 0.55. Consulting Table III we find 0.32 < 
Xx < 0.64 in contrast with typical Ax ~ 1 predictions when no SUSY corrections 
are considered. A low Ax implies that the RGBs which govern the Yukawa evolution 
are dominated by the gauge terms, thus the effect of decoupling the right-handed tau 
neutrino is small. 




10' 10^ lO'^ 10^ 10'' 10" 10" lOi'^ 10^^ 10" 
Q / GeV 



FIG. 10. Running of the third family t — b — T — Vr Yukawa coupUngs between the 

unification scale Q = Mx and the effective supcrsymmetry breaking scale Q = Ms- 
The decoupling of the tau right-handed neutrino at Q = Mj, mainly affects the 
running of (dashed line.) 

In Table IV we list the predicted values for the heavy right-handed neutrino mass 
scale Ml, which in our model is fixed by the requirement that the light left-handed 
neutrino tau has a mass of — 0.05 eV. We see that is in the range 2 — 5 x 10^^ 
GeV. 
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TABLE IV. 



Case 


A 


B 


C 


D 


E 


F 


G 


H I 




2.97 


3.35 


1.98 


2.18 


4.16 


4.77 


2.96 


3.35 4.09 



TABLE IV. Predicted values for the mass of the heavy right-handed neutrino M^, 
(given in GeV) required to generate a Ught left-handed neutrino mass m^^ = 0.05 eV. 
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FIG. 11. Running of the tau neutrino mass from the unification scale to low energy 
for Cases A,C,E,G in Table 1. 

In Fig. 11 the evolution of rav^ from Mx to Ms is shown for cases A,C,E,G. For 
each hne a different Mj,{Q — Mx) was chosen in order to assure the same rrii, — 0.05 
eV at low energy. For the purpose of illustration we have in Fig. 12 (and only in this 
plot) relaxed the constraint upon M^ and set it to vary at ten log intervals from Mx 
to lO^^Mx. The effect of changing M,^ is shown in the Mf-Xx plane. Comparing this 
plot with the tables throughout this article we can conclude that the uncertainties 
attached to cts, Mi, and Ms are far more important than the ones which might affect 
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the determination of My. The predictions for Mt and \x increase for decreasing Mi,. 
The variation in Mt is quite small (typically less than 1 %). On the other hand 
variations in Ax are larger and of about 3-4 %. 
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FIG. 12. Variations in the predictions for the top mass Mt and in the unified Yukawa 
couphng \x induced by changing M^, in the range (1— 10~^) x Mx- The dark circle 
refers to the choice of My which for Case I gives — 0.05 eV. 



VI. YUKAWA UNIFICATION IN THE 

5f/(4) SU{2)l SU{2)r model 

In this section we will extend the analysis of the top quark mass, so far considered 
in the context of the MSSM with a right-handed tau neutrino, by assuming that 
above the unification scale, at which the gauge couplings meet, the gauge group is 



S'f/(4) (g) SU{2)l ® SU{2)r |T2|. We will proceed in three steps. Firstly we review 
the model constrained by universal soft SUSY breaking masses (USBM). Secondly, 
the model is considered under non-USBM conditions and thirdly we will allow for 
the presence of D-terms which generally arise from the reduction of the gauge group 
rank ||2^. Finally, we conclude by showing that the 422 model with non- universal 
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family dependent A-terms, is compatible with third family Yukawa unification (which 
requires a negative sign for /i) and with the experimental branch ratio for the 6 — 57 
and r — > /i7 decays. 



A. The model 



Here we briefly summarize the parts of the model that are relevant for our analysis. 
For a more complete discussion see Ref. The SM fermions, together with the 



right-handed tau neutrino, are conveniently accommodated in the following F'^ = 
(4, 1, 2) and F = (4, 2, 1) representations : 

„c _ ( d'^ ^'^\ _ / u u u h'\ , . 

^~U' u^J^ d d ejj, ^ ^ 

The MSSM Higgs bosons fields are contained in h = (1, 2, 2) : 

whereas the heavy Higgs bosons H = (4, 1, 2) and H = (4, 1, 2) are denoted by : 

jj _ f Hd Hd Hd He\ jj _ f Hd Hd Hd He\ ^-^^^^ 
\Hu Hu Hu J \Hu Hu Hu Hy ) 

In addition to the Higgs fields in Eqs. (|1^) and (|1^ the model also involves an SU (4) 
sextet field D = (6, 1, 1) = {D3, Dl). [] 

The superpotential of the minimal 422 model is : 

W = F'XXABFBh + XhShh + 

\sS{HH - Ml) + XhHHD + XnHHD + FIX'j.^F-'b ^ (18) 

ivi p 



^In fact, since we wish to keep the gauge couphngs unified above Mx we also trivially include 
another pair of heavy Higgs boson fields in the (4, 2, 1) and (4, 2, 1) representations and six more 
replicas of D's. This is sufficient to assure that, above the unification scale, the one loop RGEs beta 
functions of the gauge couplings are equal. 

^Note that FFD and F'^F'^D terms, associated with baryon number violating processes, can be 
forbidden by imposing a global R-symmetry. |47| 
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where S denotes a gauge singlet superfield, the A's are real dimensionless parameters 
and Mh ~ Mx- Additionally, the Planck mass is denoted by Mp ~ 2.4 x 10^^ GeV. As 
a result of the superpotential terms involving the singlet S, the Higgs fields develop 
VEVs {H ) = {H^) ~ Mx and {H) = (H^) ~ Mx which lead to the symmetry 
breaking : 

SU{A) ® SU{2)l ® SU{2)r SU{3), ® SU{2)l ® f/(l)y. (19) 
The singlet S itself also naturally develops a small VEV of the order of the SUSY 



breaking scale so that the XhS term in Eq. (0) gives an effective fi parameter 
of the correct order of magnitude. Under Eq. (^) the Higgs field h in Eq. ([T6|) splits 
into the familiar MSSM doublets and whose neutral components subsequently 
develop weak scale VEVs = V2 and {h'^) = Vi with tan/? = V2/V1. The neutrino 
field u'^ acquires a large mass M^, ~ X'{HH)/Mp through the non-renormalizable 
term in W which, together with the Dirac u'^ - u interaction (proportional to X{h^)), 
gives rise to a 2 x 2 matrix that generates, via see-saw mechanism, a suppressed 
mass for the left-handed neutrino state The D field does not develop a VEV 

but the terms HHD and HHD combine the colour triplet parts oi H, H and D into 
acceptable GUT scale mass terms p6 |. 



In addition to the terms generated by the superpotential of Eq. ([T8|) the lagrangian 
of the 422 model also includes the corresponding trilinear soft terms F'^AFh + X^Shh, 
[| masses for the S'f/(4), SU{2)l, SU{2)r gauginos M4, M2L, M2R and explicit soft 
masses for the scalar fields rhp\F\'^ + rhpc\F'^\'^ + h\^. 



Finally we remind that the symmetry breaking in Eq (|T9|) leads to specific relations 

between the S'[/(4), SU{2)2r and U{1)y gauge couplings and gaugino masses at Mx 
''Often re-parametrised hy A = AX and m\ — \h{S). 
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given by |^ : 



5 _ 2 ^ 3 

ai a2R 



5Mi 
tti 



2M4 ^ 3M2R 



a4 



Ol2R 



(20) 



B. The universal model 



In this section we briefly review the 422 model with USBM. The main motivation 
for USBM is the smallness of flavour changing neutral currents (FCNC) and it's 
simplicity — few input parameters are needed to specify a otherwise complex and 
largely unconstrained parameter space. 

The model is very similar to the one presented in the previous section except that 
we have imposed the high energy boundary conditions at the Planck scale instead of 
at Mxi thus we allowed for the running of the parameters between these two scales. 
The relevant input parameters were sign < 0, = and : 

< Ml/2 < 2 TeV < mo < M1/2 (21) 

where M1/2 = M4 = M2L = M2R is the common gaugino mass and = fh^ = 
mjpc = ml is the universal scalar mass (a^ = 0.120, Mb = 4.8 GeV). The results are 
presented in scattered plots in Fig. 13 and Fig. 14. 

In Fig. 13 we plot ttiq against Mi/2- The black (white) circles indicate points 
that do (not) satisfy EWSB, i.e. ml — m\ < M|. We observe that mcreasmg mo is 
disfavoured by EWSB whereas increasing M1/2 makes EWSB easier to occur. The in- 
terplay between these two dependencies dictates that only the region with M1/2 > mo, 
with a threshold for M1/2 ~ 500 — 600 GeV, is allowed. The large gaugino mass 
(required in models with USBM and large tan (3) has important implications, it leads 

to a large negative m2 parameter at low energy. Thus the EWSB condition in Eq. (^) 

®See appendix D for a comprehensive derivation of these relations. 
^This condition can be obtained from Eqs. (H) and (|^) in the hmit of large tan/?. 
'^'^Note that this value corresponds to the gaugino mass at the Planck scale which in our model 
decreases between Mp and Mx- At the GUT scale we obtained Afi/2 > 400 GeV. 
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can only be satisfied if |m|| — /^^ ~ ^1/2 with |m|| ~ ^ ^1/2 — fine tuning 



(9(500) is inevitable [pO| , |49| . Furthermore large values for /x, which is correlated with 
Ml/2, increase the SUSY correction to the bottom mass in Eq. and Eq. pushing 
the prediction for the top mass below the experimental lower bound. 
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FIG. 13. Common scalar mass toq against unified gaugino mass fixed at the 

Planck scale. The white circles denote points that fail to satisfy the electro-weak 
symmetry breaking conditions. 

In Fig. 14 we show the top mass prediction Mt against Mi/2- We can see that 
not only is Mt smaller than 170 GeV (even for M1/2 as large as 2 TeV) but also that 
for fixed M1/2 its dependence with mo is small. It is relevant to note that, in the 
attractive scenario of lowest M1/2, corresponding to lighter sparticles and smallest 
fine tuning, the largest top mass prediction is unacceptable. Indeed, we found, by 
pushing the stability of the Higgs potential to the extreme, that for a low gaugino 
mass Ml/2 ~ 600 GeV □ we obtain Mt ~ 157 GeV (tan/5= 43). Nevertheless, the 
fact that in Fig. 14 the white circles are above of the blacks suggests that, if the 
scalar masses rhp^ rhpc, fhh are allowed to increase (and/or split) and be bigger than 
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Corresponding to A/1/2 (Mx) = 485 GeV and gluino masses mg{Mz) — 1145 GeV. 
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Ml/2, then Mt may increase. Thus, the main problem is how to concihate EWSB 
with scalar masses bigger than gaugino masses. 
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FIG. 14. Top mass prediction Mt against the common gaugino mass M1/2 fixed at 
the Planck scale {as = 0.120.) 



C. The non- universal model 



Over the past decade the numerous studies of models with USBM M, O, [T3, pO 



5T| , [5^ , p3[ at the GUT scale showed that they often face recurrent difficulties with 
low energy phenomenology. The unsatisfactory results have stimulated the interest 
in non-universal models |jl9|, |2l|, ^ |55|, ^ which are well motivated from 



a theoretical point of view. As was emphasized in ||55| , |56[| , even when universal- 
ity is imposed at the Planck scale, radiative corrections between Mp and Mx and 
heavy threshold effects lead to non-universal parameters at the GUT scale. More- 
over, the most general SUGRA models (with non-canonical kinetic terms) in which 
supersymmetry is broken in a hidden sector, and/or superstring theories in which 
supersymmetry is broken by the F component of the moduli fields with different 
weights, show that non-universality can be generated at the Planck or string scale. 
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In this section we consider the 422 model with non-universal boundary conditions 
at Mp. The independent input parameters were fh?p, fh?pc (proportional to the unit 
matrix), and M4 = M2R 7^ M2l- Q These parameters were made to vary at 
random in the ranges : 

200 GeV < M4 = M2R < 2 TeV (22) 

\Mi < M2L < 2M4 (23) 

200 GeV <mp^ rhpc ^ rhh < M4. (24) 

The results are presented in Figs. 15-21 which main purpose was not to exhaustively 
scan all the parameter space but to be a guide of the configuration of the parameter 
space which most effectively increased Mt 

In Fig. 15 we plot the top mass against M4. Comparing with Fig. 14 we observe 
that while the average increase of Mt with M4 is similar, the strict correlation present 
in the universal case is replaced, in the non-universal model, by a dispersed region 
of enhanced/suppressed Mt predictions. This plot also illustrates that one can ex- 
pect Mt ~ 170 GeV (for M4 > 1000 GeV) and shows that EWSB can occur under 
more relaxed conditions, which is obvious from the way the white and black circles 
distribute evenly. 

In Fig. 16 M2L is plotted against M4. We see that M2L > M4 is disfavoured by 
the condition on the Higgs potential to be bounded. The majority of black circles 
is concentrated in the M4 > M2L region with M4 > 600 GeV. This preference is 
illustrated by the left-right asymmetry around the x-axis of Fig. 17 which displays 

Mt against M4 - Msl- 

^^Note that, for the sake of simphcity, we have taken the S'C/(4) and SU {2)b, gaugino masses equal 
at Mp. Since their one- loop beta functions are identical, and the gauge couplings are roughly unified 
above Mx we also have M4{Mx) = M2r{Mx). 

^^A similar graph is obtained when Mt is plotted against M2L thus we spared from including it 
here. 
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FIG. 15. Top mass prediction Mt against the SU{4:) gaugino mass M4 fixed at 
the Planck scale in the non-universal model. The thick central line indicates the 
prediction for Mt in the universal model of Fig. 14 (a^ = 0.120.) 




FIG. 16. The SU{2)l gaugino mass M2L against the SU{A) gaugino mass M4. 
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FIG. 17. Top mass prediction Mt against the difference between the SU{4) and 

SU{2)l gaugino masses M4 — M2l- The stronger concentration of black circles on 
the right-half part of this figure indicates that EWSB favours M4 to be bigger than 

M2L- 
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FIG. 18. The top mass prediction Mi against the left-handed SUSY scalar mass 
mp- The lower mp the bigger the spread in Mt [ols = 0.120.) 
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FIG. 19. Top mass prediction Mt against the difference between the left/right SUSY 

scalar masses tjif — m^c . This symmetric plot shows that the top mass does not 
favour an hierarchy between rfiF and mF" (ccg = 0.120.) 
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FIG. 20. Top mass prediction Mt against the soft Higgs mass fhh- We can see that 
a large top mass is favoured by a large Higgs parameter rhh (ccs = 0.120.) 
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FIG. 21. Top mass prediction Mf against the ratio between the soft Higgs mass rhh 
and the average of the SUSY scalar masses ^{mp +mF<^)- The bigger this ratio the 
bigger Mt is hkely to be. 

The Fig. 18 shows that, generally, as itif decreases Mt predictions are allowed to 
increase. A similar dependence of Mt on mi?c was found. Indeed, when Mt is plotted 
against the difference Arhp = rhp — rhpc such as in Fig. 19 a symmetric left-right 
graph emerges around the x-axis. On the other hand, from Fig. 20 we can conclude 
that Mt increases with increasing soft Higgs mass rhh- These figures suggest that 
the top mass may be increased by increasing the splitting between the sfermions and 
Higgs soft mass. We therefore find interesting to plot in Fig. 21 against the ratio 
"nih/ (t^f) where {niF) — {'mF + 'mFc)/2 is the average sfermion mass. Our suspictions 
are confirmed. 

After combining the top mass dependencies on the input parameters suggested 
by the previous figures with numerous case-by-case analyses of "outputs" from our 
numerical model we arrived at the following conclusions. The top mass prediction is 
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strongly dependent on M4. Q Once M4 is fixed, decreasing fhp and Thp': increases 
both Mt and the Higgs potential stability parameter (that is positive for a "bounded- 
from-below" Higgs potential) : Q 

= + fil - 2\ml\ > (25) 

However rfippc cannot be arbitrarily small since if rfip^pc are very small the sleptons 
became too light. Increasing M2L and rhh increases the top mass, however they 
also affect S. We found that the most efficient mechanism to increase Mt relied on 
decreasing M2L (which decreased Mt moderately but increased S substantially) and 
increasing rhh (which increased Mt significantly and decreased 5* moderately). For the 
sake of illustration, we found that, by pushing the EWSB conditions to the extreme 
{S ~ 0, I.e. tuning mg), it was possible to get Mt ~ 175 GeV only if M^^Mp) > 800 
GeV. In this case we got M2L = 500 GeV, fhp = 350 GeV, rfipc = 450 GeV and 
rhh = 1000 GeV. Q Naturally, for a larger M4 mass, S increases and an acceptable 
top mass is easier to obtain with less tuning in 777,3. 

The main conclusion is that in the 422 model with non-universal soft masses a 
top mass of around Mt ~ 175 GeV is only possible to obtain in the context of a large 
gluino mass mg{Mz) ~ 1520 GeV and large tan/5 ~ 50, by implicitly tuning the 
EWSB conditions, and by choosing squark/slepton masses considerably smaller than 
the soft Higgs mass rhp,rhpc < rhh. 



-'^''The reason is because, at the unification scale, the gluino masses are set to be Ms{Mx) = 
M4^(Mx)- Thus the sensitiviness of the results with M4 is in fact a sensitiviness to the masses of 
the coloured sparticles. 

^^In the large tan/3 limit S"^ is given by S*^ ~ to^ — TOj — A/| ~ to§/ tan/3. 

i^The corresponding values at Mx are M4 = 653 GeV, M2L = 408 GeV, mF = 456 GeV, 
TOF<= = 567 GeV, and fhh = 972 GeV. 
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D. D-term contributions 



The symmetry breaking of the Pati-Salam to the SM gauge group given in Eq. ( |19|) 
reduces the group rank from five to four. Although the GUT symmetry is broken 
at a very high energy it nevertheless has important consequences to low energy TeV 
phenomenology via D-term contributions to the scalar masses ^ |5^. In the 422 
model the GUT boundary conditions for the scalar masses are : Q 



ml 


= ml + glD^ 




ml. 


= m% - {gl - 2gl^] 




m% 


= m% - {gl + 2gl^] 




mf 


= ml-3glD^ 




m% 


= Th% + {3gl-2gj^)D'' 


ml 


= ml. + {3gl + 2glj,)D^ 


ml 


= ml-2gl^D^ 




ml 


= ml + 2gl^D^ 





(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



The D-term corrections, If D is sufficiently large, are important to consider be- 
cause, firstly they leave an imprint in the scalar masses of the charges carried by 
the broken GUT generator (these charges determine the coefficients of the g"^ terms 
above), therefore the analysis of the sparticle spectra |Q might reveal the nature of 
the GUT symmetry breaking pattern; Secondly, they split the soft Higgs masses by 



rrin — m 



-AgxD"^, which for positive D^, makes radiative EWSB much easier to 



occur. Indeed, we found that once D > 150 — 200 GeV then EWSB no longer requires 

the large gluino/squark masses characteristic of models with USBM. 

^^See appendix E for a detailed derivation of the Z?-terms in the 422 modeL Note that in the hmit 
of unified gauge couphngs the D-terms give identical corrections to the fields in the 16 dimensional 
representation of 5*0 (10). 
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FIG. 22. Correlation between the sfermion masses •fhp = niF^ and the Higgs soft 
mass rhh required to obtain a top mass around 175 GeV {as = 0.120.) Three il- 
lustrative choices for the value of the D-term are taken. The effect of experimental 
uncertainties in the top and bottom quark masses is indicated by the toned ar- 
eas. The box around the D = 200 GeV label indicates the scanned region used to 
generate Figs. 23-31 (see main text for details.) 

In Fig. 22 we summarize how an acceptable top mass prediction can be achieved 
in the 422 model in the 5*0(10) limit of universal gaugino masses M1/2 and universal 
squark and slepton mass parameters itif = rhpc defined at the Planck scale, but 
with D-teim. corrections arising at the symmetry breaking scale. Q The three solid 
lines establish the correlation between the sfermion masses rhp = rhpc (plotted in the 
y-axis) and the Higgs soft mass rfih (plotted in the x-axis) required to obtain a top 
mass Mt = 175 GeV, for D = 200, 300, 400 GeV and M4 = M2R = M2L = 300 GeV. 
For the D = 200 GeV line, EWSB failed to occur beyond the point where the line is 
discontinued in her way upward. In all cases the lines were cut at the bottom edge 

at points where the left-handed stau mass became smaller than 100 GeV. 

-"^^ Although we refer to this as the S'O(IO) limit, in fact the two theories differ in the region between 
the GUT scale and the Planck scale. 

^^The D — 300,400 GeV lines continue beyond the border of the figure. 



35 



The general behaviour in Fig. 22 with increasing D-teim values — lines shifting 
to the top-right corner — results from the decrease of m| ~ rh'j^ — 3gxD^ in Eq. ( p9D 
which demands an increasing parameter. The broader darker areas around the 
solid lines correspond to the uncertainty in the top mass 175 ± 5 GeV, while the 
thinner patches result from varying the pole bottom mass Mb in the range 4.8-5.1 
GeV. We also analysed how the gaugino masses affect the solid lines. For example, 
for D = 200 GeV, taking M1/2 = 400 GeV generates the dashed line. If M2R or M2L 
are taken different from M4 the prediction for the top mass is very similar to the one 
with universal gaugino masses with M1/2 = M4. 

In summary. Fig. 22 shows that if D> 200 GeV then a successful prediction for 
Mf can be achieved, even for small gluino/squark masses — no fine tuning in EWSB 
conditions — as long as the soft Higgs mass rhh is bigger than the sfermion masses 
rhp, rfiFc As an example, setting D = 200 GeV and M1/2 = 300 GeV we found that 
Mf ~ 175 GeV could be obtained by fixing rfiF = fhpc = 500 GeV and rhh = 750 
GeV. This input lead to the following low energy predictions : tan /5=52, Ax = 0.70 
a gluino mass rrig = 618 GeV, a lightest neutralino mass m^o = 88 GeV (bino like, 
but with substantial higgsino component), a lightest chargino mass m^- = 125 GeV 
(higgsino like) and masses for the lightest sfermions mc,^ = 179 GeV, m^^ = 189 GeV 
and = 377 GeV. The lightest CP-even Higgs mass, computed using the one-loop 
expressions of Ref. (that include the stop/sbottom corrections only) was found 
to be rrih = 114 GeV (note that this value should be read with some caution, we 
estimate an error of about 10 GeV in rrih.) 



2° The corresponding values at Mx are M1/2 = 247 GeV, rfiF.F- = 474 GeV and rhh = 691 GeV. 
After the D-term corrections the scalar masses are : rhg — rhuc = TOgc ~ 496 GeV, m; = to^c = 404 
GeV, m^c = 572 GeV (for the third family) and = 660 GeV, mi = 720 GeV. 
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In the last part of this section we present a series of figures that show in detail how 
an acceptable top mass prediction can be achieved and how it is correlated with the 
sparticle spectrum. For the sake of illustration we take, as in the previous sections 
tts = 0.120, Mfc = 4.8 GeV and the following input at the Planck scale : 

Ml/2 = 300 GeV D = 200 GeV 

(34) 

rhh = 700 - 1100 GeV rhp = rhpc = 500, 550, ..,900 GeV 

These values correspond to a scan of the rhh-rhp^pc parameter space that is indicated 

in Fig. 22 with a box. 

In Fig. 23 we plot the pole top mass prediction Mt against rhh for several choices 
of the sfermion masses rfip^pc = 500, 550, .., 900 GeV. We observe that Mt increases 
with increasing rhh and decreasing rhp^pc The reason for such dependence is directly 
related with the value of the Higgs mixing parameter /i. 

In Fig. 24 we show how the top mass prediction Mt correlates with fi. In this plot 
each line corresponds to a fixed fhp^pc mass (labeled by the number beside it) and 
along it ihh is varying in the range indicated by Eq. ( p4D . We see that decreasing m p^p<^ 
shifts the lines to the right of the graph thus making /i less negative. Furthermore, 
as rhh increases, from the bottom to the top of the graph, decreases. These 
dependencies lead to a small value for /i (at the top-right corner of this graph) which 
in turn lead to small SUSY bottom corrections that raise the top mass prediction to 
an acceptable value. 

In Fig. 25 we show the correlation between the top mass Mt and the gluino mass 
rrig predictions. We see that, for the choice of M1/2 = 300 GeV in Eq. (p4D, if the top 
mass is in the 170-180 GeV range then the gluino mass is in the 605-620 GeV range. It 
is worth stressing that the 422 model with D-terms offers the possibility of predicting 
light gluinos. This is exciting from the experimental point of view and theoretically 
desirable since it reduces the fine-tuning in the Higgs potential parameters. 
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FIG. 23 Prediction for the top mass Mt against the soft mass for the unified Higgs 
bosons fhh in Eq. (p^. The numerical labels, for each line, indicate the value for 
the soft masses of the scalar fields rhp = rhpc used to generate the corresponding 
line {as = 0.120.) 
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FIG. 24 Prediction for the top mass Mt against the Higgs mixing parameter /i. Each 
line corresponds to a fixed choice for the soft mass of the scalar fields •mp — mF" and 
along it the soft mass for the unified Higgs bosons rhh is varying (increasing from 
the bottom to the top of the plot.) We observe that Mt increases with decreasing 
ImI (a, = 0.120.) 
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FIG. 25 Correlation between the top mass prediction and the gluino mass pre- 
diction nig. Each Une corresponds to a fixed choice for the soft mass of the scalar 
fields TTiF = rriF" (labeled with a number) and along it the soft mass for the unified 
Higgs bosons m/j is varying (increasing from the bottom to the top of the plot.) 




FIG. 26 Correlation between the top mass prediction Mt and the; lightest chargino 
mass prediction m - . Each line corresponds to a fixed choice for the soft mass of 
the scalar fields rfiF = rhF" (labeled with a number) and along it the soft mass for 
the unified Higgs bosons rhh is varying (increasing from the bottom to the top of 
the plot.) 
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In Fig. 26 we plot the top mass Mt against the lightest chargino mass m - pre- 
diction. In this model Xi is roughly the charged higgsino which is lighter than 
the charged wino because < at low energy. From this graph we read that 
80 GeV < m - < 180 GeV which should be compared with experimental bound 



from LEP2 : m'!.^ > 89 GeV 



In Fig. 27 we show the correlation between the top mass Mf and the lightest 
neutralino mass m^o predictions. We note that since /i can be comparable with of 
the bino mass, Xi can have a substantial higgsino component. From this plot we read 
that 60 GeV < m^o < 105 GeV which should be compared with the experimental 
bound : > 40 GeV 

In Fig. 28 we plot Mj against the lightest charged slepton mass m^^ prediction. It is 
interesting to observe that the D-term correction to the left-handed charged sleptons 
in Eq. (^) is negative while for the right-handed charged sleptons in Eq. (|30D it is 
positive. Thus the lightest stau is not right-handed (f2) but left-handed (fi). From 
this graph we see that the prediction for m^^ can vary significantly, roughly we find 
that only an upper bound can be imposed m^^ < 600 GeV (the experimental lower 
bound is > 81 GeV. m) 



In Fig. 29 the prediction for the top mass Mf is plotted against the lightest sneu- 
trino mass rrij^^. We note that for the choices of rhppc = 500, 550, 600, 650, 700 GeV 
the Ur mass is driven to zero as rhh increases from the bottom to the top of the graph 
(thus, it is possible that could be the lightest SUSY particle (LSP).) Comparing 
this figure with Fig. 28 we see that the prediction for the fi and 9^ masses are similar. 
We find that the predicted upper bound for m/>^ < 650 GeV is compatible with the 
experimental lower bound ml^^ > 43 GeV ^2 . 
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FIG. 27 Correlation between the top mass prediction Mt and the hghtest neutrahno 

mass prediction rn^o . Each hnc corresponds to a fixed choice for the soft mass of 
the scalar fields rriF = rriFo (labeled with a number) and along it the soft mass for 
the unified Higgs bosons rhh is varying (increasing from the bottom to the top of 
the plot.) 




FIG. 28 Correlation between the top mass prediction Mt and the lightest charged 
slepton mass prediction mf-^ . Each line corresponds to a fixed choice for the soft 
mass of the scalar fields mp = mp<: (labeled with a number) and along it the soft 
mass for the unified Higgs bosons fhh is varying (increasing from the bottom to the 
top of the plot.) 
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FIG. 29 Correlation between the top mass prediction Mf and the hghtest sneutrino 

mass prediction m.^,^. Each line corresponds to a fixed choice for the soft mass of 
the scalar fields mp = rriFc (labeled with a number) and along it the soft mass for 
the unified Higgs bosons rhh is varying (increasing from the bottom to the top of 
the plot.) 




FIG. 30 Correlation between the top mass prediction Mt and the lightest sbottom 
mass prediction mg^. Each line corresponds to a fixed choice for the soft mass of 
the scalar fields rhp = fiiF^ (labeled with a number) and along it the soft mass for 
the unified Higgs bosons rhh is varying (increasing from the bottom to the top of 
the plot.) 
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In Fig. 30 the top mass prediction Mt is plotted against the hghtest squark mass, 
that is, the right-handed sbottom 62- The reason why the right-handed down-type 
squarks are hghter than the other squarks is because of the negative D-term correction 
in Eq. (PB]). The prediction for m^^ is well above the experimental lower bound 



mf^ > 75 GeV 
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FIG. 31 Correlation between the top mass prediction Mt and the CP-odd Higgs 
boson mass prediction niA- Each Hne corresponds to a fixed choice for the soft mass 
of the scalar fields ifip = ifiF" (labeled with a number) and along it the soft mass 
for the unified Higgs bosons ihh is varying (increasing from the bottom to the top 
of the plot.) 

Finally in Fig. 31 we present the correlation between the predictions for the top 
mass Mt and the mass of the CP-odd Higgs boson itia- We observe that decreases 
for increasing rhp^F'^ masses. In fact, for fhp^pc = 900 GeV we find that is driven 
to zero. This graph shows that only an upper bound for the Higgs boson mass can be 



estimated : < 220 GeV (the experimental lower bound is m^^^ > 84 GeV 



We would like to end this section by remarking that, as far as the top mass 
prediction is concerned, the results in this section are similar to the SO{10) model 
0, since we have taken the case of universal gaugino and scalar 
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mass parameters, with the D-term corrections reducing to those in SO (10). However 
we have included the neutrino Yukawa couphng in our analysis, and also the theory 
is different from 5*0(10) above the GUT scale. The reason we worked in this limit 
was to make contact with the recent work on 5*0(10), and to distinguish clearly the 
effects of D-terms from the effect of explicit non-universality. However we would like 
to emphasise that the results in Figs. 23-31 are not only valuable by themselves, but 
they also serve as a reference to which the results in section VI. C should be compared 
with. Generally, the non-universal 422 model with D-term corrections will lead to a 
set of predictions that is a combination of the results presented in sections VI. C and 
VI. D. 

E. 6 ^ S7 and r ^ /i7 
with non-universal A-terms 

It is well known that the decay 6 ^ S7 is a sensitive probe of new physics. In 
the standard model the loop diagram involving the W boson and top quark give a 
theoretical prediction for BR(6 ^ S7) = (3.28 ± 0.33) x 10"^ H □ at the next- 
to-leading order in QCD. This result turns out to be slightly larger than the official 
CLEO measurement (see first Ref. in |6^). However the recent ALEPH results 
indicate a larger branching ratio [Q. We quote : 



BR(6->S7) = (2.32 ± 0.57 ± 0.35) X 10-^ CLEO 

(35) 

BR(6^S7) = (3.11 ±0.88 ±0.72) X 10"^ ALEPH 

where the first error is statistical and the second is systematic. In view of the above 
large uncertainties we will take the conservative range : 

1.0 X 10-^ < BR(6 57) < 5.0 X 10"^ (36) 



21 This is the value of K. Chetyrkin et al. in Rcf. |2[. 

Note that an updated prcUminary value by CLEO has been reported : BR(& 57) = (3.15 ± 



0.35 ± 0.32) X 10"'' (see second Ref. in 
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In supersymmetric extensions of the SM the inclusion of additional SUSY par- 



ticles can spoil the above "agreement" The reason is because, generally, there 



is no guarantee that the sparticle contributions are small or that they conspire to 
cancel between themselves. Indeed, it is known that in two Higgs extensions of the 
SM the contribution from the charged Higgs-top quark loop H^t always interferes 
constructively with the SM one. Thus, the resulting branch ratio prediction is bigger 
and dependent on the unknown charged Higgs mass rriH-- If ttlh- is light then the 
prediction for the BR(6 — 57) is larger than the experimental upper bound. On the 
other hand, for a heavy Higgs mass rriH- ~ 1 TeV, agreement with Eq. (|35|) is still 
be possible. 

The above situation can change drastically in SUSY models with third family 
Yukawa unification. The reason is because the new contributions arising from gluino- 
sbottom (gb), neutralino-sbottom (x^b) and particularly from chargino-stop (x t) 
loops are enhanced by large large tan /5 ~ 40 — 50 factors. In contrast to the charged 
Higgs contribution, the SUSY amplitudes can either add constructively or destruc- 
tively with the H~t contribution or between themselves. In order to recover the 



original agreement with experiment a large part of the work in the literature p5 
explores the idea of suppressing the 6 ^ S7 decay by canceling the x~t against the 
H^t amplitude (typically the gluino and neutralino amplitudes are small.) Many of 
these results are obtained in the context of SUGRA/GUT models with universal soft 
masses and trilinear terms at the Planck or GUT scale with large gaugino masses 
Ml/2 ~ 1 TeV. Cancellation between x~i and H~t loops is possible because of two 
reasons, firstly the large gaugino mass leads to large stop masses thus suppressing 
the x~t amplitude making it comparable in magnitude to the H~t term. Secondly, 
the sign of fi is chosen to be positive. This means that the x~t contribution has the 
opposite sign of the H~t contribution thus allowing the cancellation to occur. 

Clearly the above strategy is very attractive because it is based on universal soft 
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parameters which render it model independent. However, it fails to address two issues. 
Firstly a large M1/2 leads to fine-tuning of the Z boson mass and secondly, a positive 
is in conflict with the successful prediction for the top mass in the context of GUT 
models with unified third family Yukawa couplings which, as we have seen, require 
negative \i. 

In this section we explore a different solution to the enhanced b ^ decay in 
SUSY models with t — b — t Yukawa unification. The idea is to choose /x to be 
negative, for the reasons stated earlier, and to suppress 6 — > 57 by allowing the 
trilinear soft A-terms to have a non-universal family dependent structure. We will 
avoid unnatural tuning of electroweak symmetry breaking by taking a low value for 
the gaugino masses, for example M1/2 ~ 300 (250) GeV at the Planck (GUT) scale. 
In this scheme we find that 6 — > 57 is dominated by the chargino-stop loop. 

The purpose of this section is to investigate the possibility of suppressing 6 — > 57 
by tuning the initial non- universal A-terms at the Planck scale such that, at low 
energy, they are cancelled by the fiavour violating signals that naturally develop 
when the parameters of our model evolve from high to low energy by the use of 
RGEs. We will conclude by checking that the introduction of non-universal A-terms 
is also compatible with the present upper bound on BR(t — > /ij). 
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b ^ s^y 



The standard model, charged Higgs, chargino and gluino diagrams contributing 
to this decay are illustrated in Fig. 32. 




FIG. 32. Diagrams responsible for the 5 ^ 37 decay considered in this article (note 
that we did not compute the neutralino/sbottom loop which is typically found to 
be small.) The W boson loop a) illustrates the standard model contribution. The 
charged Higgs loop b) is present in two doublet Higgs extensions of the SM. The 
gluino contribution c) was found to be small. On the other hand, the chargino loop 
d) dominated the b sj decay. 



In models with large tan f3 the corresponding dominant amplitudes are given by : 
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(37) 
(38) 
(39) 

(40) 



where = g'^/^An), ae = e^/(47r), M]y is the W boson mass, m^^ and m^^ 
the sup/sdown type squark mass eigenstates, V the CKM matrix, xtw 



m't/M^, 



Xth = m^/mjj-, X 



The flavour matrices V^'^ describe the mismatch between the transformations that 
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diagonalize the down-type squark and quark mass matrices (See appendix A for de- 
tails.) The functions F are defined by Fmnix) = ^Fm{x) + Fn{x) with Fm^n{x) as in 
Ref. [^. The gaugino and higgsino vertex factors G and H are given by : 

= T^AV^%A (41) 
Hijf = T^2j:{V}rU{Yu)BBVBA (42) 



B=l 



Hht = S^^^{Y,)UVr)A. (43) 



where, 



Yu = diagim^, mc,mt)/{V2Mw sin P) = X'Jg (44) 
Yd = diag{md, TUs, nib) /{\/2Mw cos (3) = X'Jg. (45) 

and S'^, diagonahze the chargino mass matrix. Finally the precise definition of 
the V^'', V^" matrices that describe the mismatch between the transformations re- 
quired to diagonalize the up-type squarks and quarks mass matrices can be found in 
appendix A. 

The branch ratio BR(6 — » 57) can be computed from [^] : 

BR(6 ^ 57) = 5.?^^^I bR(& ^ cev) (46) 
r(o — i> cez/) 

where the decay rate for 6 ^ 57 is given by : 

T{h^s^) = ^\A\mb)\\ (47) 

The full expressions for T{h ^ cez/) and A^{m}y) (the QCD corrected amplitude 
at the scale of the process (~ m^) - obtained from the total sum of amplitudes 
A^{Mw) = Asm + Ah- + A^- + Ag) can be found in Ref. 



It is interesting to note that the chargino amplitude has two distinct contributions. 
They are illustrated in Fig. 33. In Fig. 33 a) the helicity flip required in the decay is 
achieved at the higgsino vertex. Thus, along the internal squark line ti-ci, flavour 
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violation develops through a (A23)ll mass insertion. On the other hand, in 
Fig. 33 b) the two higgsino vertices require that the helicity flip must be accomplished 
through the t^-CR line via a {^'^■ilLR mass insertion. 




FIG. 33 Diagrams corresponding to the two contributions associated with the chargino 
amplitude of Eq. (p9|). In fig. a) flavour violation is introduced through a (A23)ll 
mass insertion along the t^-CL squark line. In fig. b) a {/^2z)lr mass insertion is 
introduced along the Il-cr chirality flipping squark line. 

Remembering that we aim at reducing A^- by introducing tree-level explicit 
sources of flavour violation through non-universal soft SUSY terms, the relevant ques- 
tion to address is whether it is more appropriate to relax universality of the soft mass 
terms for the fields in the F, F'^ multiplets, thereby changing (A23)ll, or to modify 
the trilinear A-terms, thereby changing (A23)lh- Foi' reasons that will became clear 
latter, when we study the r — /i7 decay, we will take the latter approach. 

Traditionally the SUSY trilinear terms are defined, at high energy, to be propor- 
tional to their associated Yukawa matrix, i.e. {Au)ab = Ao{.^u)ab etc.. Although 
this assumption leads to obvious simplifications, models inspired by string theory 
have been proposed in which the A-terms are not universal. Motivated by these 
results we parametrise the A's by : 

{Au)aB = AoXab{^u)aB i^dlAB = ^oXAsi^djAB 

(48) 

{Ae)AB = AoXab{K)aB {^u) AB = AoXab{K)aB 

where xab is a dimensionless matrix of order one that we conveniently choose to have 

^■^The (A23)ll parameter is the off-diagonal 23 entry of the left-handed up-type squark mass 
matrix in a basis where the up- type quark mass matrix is diagonal (see Ref. pH].) 
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the following structure : 



xab =11 X 




] 



(49) 



We now turn to study the phenomenological implications of Eq. ( ^9]) to the 
BR(6 — i> 57). For illustrative purposes, and keeping in mind the successful pre- 
diction for the top mass in the 422 model with Z^-terms, we use, as in the previous 
section, the same = 0.120, Mf, = 4.8 GeV and the same input at the Planck 
scale : My 2 = 300 GeV, = = 500 GeV, rhh = 750 GeV and D = 200 GeV. 
Moreover we fixed Aq = 2000 GeV and a negative sign for fi. At the GUT scale 
the Yukawa matrices were set as follows. The eigenvalues of A„, A^, Ae were fixed by 
the requirement that the fermion masses at low energy were reproduced. The CKM 
mixing angles were assumed to be given by di = 6^ + Of with 9"^ = 6f (the 6'"''^ are 
the angles that parameterise the rotation matrices that transform the left-handed 
up/down quarks into their physical mass basis.) The angles associated with the ro- 
tation required to transform the left handed charged leptons into their physical mass 
basis were given by 6^ = Of. Finally Xy was set equal to A^. For completeness we list 

in Table V the values of the Yukawa matrices at the SUSY scale Q = Ms = 430 GeV. 
24^0 (Mx) = 570 GeV 
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TABLE V. 
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539 X 10- 


-5 




\ 
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-1 


141 


X 
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TABLE V. Values of the Yukawa matrices at the effective SUSY scale Q = Ms = 
430 GeV. The correct quark and charged lepton masses as well as the CKM matrix 
are obtained after "running-down" these matrices firstly, from Q = Ms to Q = Mz, 
using the SM RGBs, and secondly from Q = Mz to Q = I GeV, using SU{3)c x 
U{l)em RGEs. (See section II. B for more details.) 
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FIG. 34 Individual values for the amplitudes contributing to the b s"/ decay 
against the trilinear parameter x of Eq. (^) . One observes that the chargino am- 
plitude A^- is sensitive to x and bigger than Asm, A^- and Ag for almost all the 
X values. Furthermore, the sign of A-^- changes as x increases from 0.0 to 2.0. 



In Fig. 34 we plot the SM, charged Higgs, chargino and gluino amphtudes in 
Eqs. (|37|)-(^0D against x. For x ~ we approximately recover the usual universal 
y4-term model. This is because the first and second family of the trilinear terms 
are small (due to small Yukawa entries) and the third family A-terms at low energy 
are not too sensitive to their initial value at high energy. We see that because /i is 
negative, the chargino amplitude is negative and interferes constructively with the 
other amplitudes. However, as x increases, the chargino amplitude becomes less 
negative and the magnitude of l^^^--] steadily decreases. At some point, around 
X ~ 1.1 the chargino amplitude vanishes. Beyond x ~ 1.1, A^- is positive, thus it 
interferes destructively with the SM and Higgs amplitudes. Clearly, we see that by 
tuning X it is easy to find a region where 6 — 57 is suppressed. 
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FIG. 35 Branch ratio for the b — > 57 decay against the trihnear parameter x in 
Eq. (^9|). The dashed hne A corresponds to the resuh one obtains when Aq = 0. 
The sohd Hne B was computed using the SM, charged Higgs, gluino and higgsino 
contributions. For the hne C we used the SM and Higgs contributions only. The 
hne D was obtained using the SM amphtude alone. The allowed experimental range 
is indicated by the area between the two dotted lines. We observe that for the line 
B an acceptable prediction is obtained for x ^ 1.60. 



In Fig. 35 we show the BR(6 —* S7) against x. The area between the dashed hnes 
indicates the experimental allowed range in Eq. (|36D. The dashed line A corresponds 
to the result one obtains when Aq = 0. The solid lines D and C are plotted for 



illustrative purposes; they correspond to the BR(6 —>■ 57) that is obtained when only 
the SM amplitude, and the SM plus the charged Higgs amplitudes are considered 



respectively. The solid line B indicates the BR(6 57) prediction when all the 
amplitudes in Eqs. (^)-(0) are considered. We observe that, due to a light charged 
Higgs mass mjj ~ 130 GeV, the value of the BR(6 — > 57) when only the SM+H~t 



loops are considered lies above the experimental range. On the other hand, the 



BR(6 s'-f) when the chargino contribution is included (line B) starts at x ~ very 
large but is driven into compatibility with experiment at around x ~ 1.60 ± 0.15. 
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In the previous section we showed that the parameterization of the trihnear terms 
by Eq. ( ^9] ) lead to a suppressed BR(6 57) compatible with the experimental 
data. However, one should be careful about other possible implications derived from 
this new source of flavour violation. Clearly Eq. (|49|) also introduces lepton flavour 
violating signals through the A^^u matrices in Eq. (^8|). The purpose of this section 
is to check that the decay r ^ fi'y (for the successful value of x fixed by 6 ^ 57) is 
not in conflict with the experimental upper bound on BR(r —>■ fi'j) < 3.0 x 10^^ . 

The effective lagrangian for the t ^ ^'j decay is : 

Cr-.^^ = l^^ip - q) [A^Pr + A^Pj) a^^Mp)Fai3 (50) 

where Prl are projection operators and F^/s the electromagnetic field tensor. The 
branch ratio can be computed using : 

BR(r-./.7) = ^(|A^P + |A^P) (51) 
where A^ = A^. + A^o and A^ = A'^o- Numerically we found that the decay is 

A. A X 

dominated by the chargino-sneutrino loop diagram illustrated in Fig. 36. 




FIG. 36 Chargino diagram involved in the t fij decay (the neutrahno diagram 
was found to be subdominant.) In this process lepton flavour violation develops 
through the (A23)ll mass insertion along the sneutrino line. 



The corresponding amplitude is given by : 
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(52) 
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where is the physical mass of the sneutrinos, -- = m? ^/rn^-, and J a dimen- 
sionless function. |^ The gaugino and higgsino vertex factors are given by : 

Git = T^iiVllU 

Hit = S^,l{Y^)AA{VrL^)Aa (53) 

where = diag(me, m^, m^)/(-\/2M(4/ cos/?) = X'^^/g. Finally the sneutrino/charged 
lepton flavour matrix is defined by Vj^^ = SI^T'^'^ where S^^^ diagonalizes the sneutrino 
mass matrix and is the matrix that rotates the left-handed charged leptons into 
their physical mass eigenstates (see appendix A for details.) 

We see from Fig. 36 that the r decay develops through the sneutrino t'J-^'^ scalar 
line via the (A23)ll mass insertion. Comparing the bottom with the tau decay, 
we note that while in the former decay the presence of right-handed up-type squarks 
allows the diagram in Fig. 33 b) to exist, in the latter decay, the corresponding lepton 
analogue of such diagram does not exist since the effective low energy theory does 
not include right-handed sneutrinos {vr}- 

We can now justify why we preferred to consider non- universal A-terms rather 

than non-universal soft SUSY masses. If we had chosen to introduce non-universal soft 

squark masses then unification of quarks and leptons in the F and F'^ multiplets would 

demand similar non-universal slepton masses. In this scenario it would be difficult 

to simultaneously suppress 6 ^ s + 7 and keep BR(r fij) below the experimental 

bound. The reason is because both decays can proceed via a (A23)ll niass insertion. 

In contrast, when non- universal A-terms are introduced, we can control 6 ^ 57 via 

the (A23)l_r insertion, while leaving the prediction for BR(r — > fi'-f) approximately 

unchanged because the leptonic decay proceeds via the (A23)ll insertion. 
25Explicitly, Jix) = [a;2 - 4a; + 3 + 2 ln(a;)]/[2(2; - l)^]. 
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FIG. 37 Values for the chargino and neutralino amplitudes A^l,^ contributing 
to the T ^ fij decay against the trilinear parameter x in Eq. . We observe that 
A^o^ are small compared with A^. Moreover, the magnitude of A^ slowly decreases 
as X increases. 
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0.0 0.5 1.0 1.5 



2.0 



FIG. 38 Branch ratio for the t fij decay against the trilinear parameter x in 
Eq. (|9|). The dashed line A corresponds to the result when Aq = 0. We observe 
that the solid line B is driven to smaller values as x increases. In fact, for large 
X ^ 2.0, the predicted branch ratio is roughly one order of magnitude smaller than 
the present experimental bound (here indicated by the dotted line.) 
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In Fig. 37 and Fig. 38 we show how the amphtudes and the branch ratio of the 
r — >• /X7 decay depend on the trihnear parameter x in Eq. P9|). We see that, for 
X ~ 0, the prediction for the tau decay is close to the experimental limit BR(r 
fi'-f) < 3.0 X 10~^. However as x increases the magnitude of the chargino amplitude 
l^-^-l slowly decreases thus driving the prediction for BR(r fi'-f) to smaller values. 
For X ~ 1.60 the predicted BR(r — > ^^7) is about one order of magnitude below the 
experimental upper bound. 

In summary, we have verified that by introducing non-universal A-terms one can 
simultaneously satisfy the present experimental constraints on 6 ^ 37 and t fij 
even when is negative. 

Finally we point that the success in suppressing the b ^ S'j decay is not without 
some tuning. Indeed we found that the suppression works well only for a restricted 
range ofx~ 1.60 ±0.15. Thus taking as a measure of tuning x the ratio 6 = 6x/x 
we find 5 ~ 19 %. Nevertheless, it is unclear if it is more natural to expect b —>■ sj 
to be suppressed by flavour physics, implying tuning in the flavour parameters, or 
by suppressing the charged Higgs and chargino amplitudes by demanding large mn- 
and stop masses, implying tuning in setting Mz- 
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VII. CONCLUSION 



We have studied Yukawa unification, including the effects of a physical neutrino 
mass consistent with the Superkamiokande observations. We began our study by 
reviewing the usual mSUGRA scenario with universal soft mass parameters, but in- 
cluding the effect of the neutrino Yukawa coupling. Assuming hierarchical neutrino 
masses, and mixing angles arising mainly from the neutrino sector (so that the charged 
lepton Yukawa matrix has no large off-diagonal entries) we saw that the usual pre- 
dictions are not much affected. Q For example the usual result that positive /i is not 
allowed, and negative /i leads to top quark masses which are too small is still valid. 

We then analysed a string/Z^-brane inspired SU{A)®SU{2)L^SU{2)fi model since 
this allows the most general non-universal scalar and gaugino masses, and is therefore 
the perfect laboratory for studying the effects of non-universality. We explored the 
sensitivity of the predictions to variations in all these non-universal soft masses, and 
showed that the usual results can change considerably as a function of the degree 
of non-universality of the soft parameters. We then switched on the usual D-teim 
contributions which arise in 50(10), and which by themselves are enough to allow 
successful electroweak symmetry breaking, and permit small corrections to the b- 
quark mass which implies an acceptable large top quark mass. We studied the effect of 
the D-terms in the SU{4:)®SU (2) l^SU (2)/j model, for the sake of clarity working in 
the (approximate) 5*0(10) limit of the model, in order to distinguish clearly the effect 
of the D-terms from that of explicit non-universality. Including D-terms as the only 
source of non-universality, we studied the sparticle spectrum in some detail, although 
these results should be considered in conjunction with the effect of the explicit non- 
universal scalar masses and non-universal gaugino masses considered previously, since 

in a realistic 5f/(4) (g) SU{2)l ® SU{2)ji model both effects will simulaneously be 

^^Large off-diagonal entries in the neutrino Yukawa matrix will not affect our results very much 
since the right-handed neutrinos decouple from the right-hand sides of the RGEs at high energy. 
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present. 

We then turned to rare decays such as 6 — > 57 and r — > //7 which severely 
constrain Yukawa unification. Assuming that the only source of non-universality 
is due to D-tcrms, wc showed how non-universal trilinear parameters can lead to 
cancellations of the important effects and so provide information about the family- 
dependent supersymmetry breaking soft Lagrangian. Again, the effect of more general 
non-universal scalar masses and non-universal gaugino masses associated with the 
SU {4)^SU (2) L®SU {2)r model will modify the results, but the main message remains 
clear: family-dependent non-universality will play an important role in 6 — > 57, and 
such effects should be considered in conjunction with the effects of family-independent 
non-universality. We have shown that successful Yukawa unification can be achieved 
if both types of non-universality are simultaneously present. 

In summary Yukawa unification is well motivated from both a theoretical point 
of view, and from the point of view of predicting large tan (3 which helps to raise the 
Higgs mass. We have found that Yukawa unification is perfectly viable providing the 
soft masses are non-universal in both family-independent and family-dependent ways, 
and we have explored the specific correlations between the soft parameters required 
in order to satisfy all the constraints simultaneously. In our view the sensitivity of 
Yukawa unification to soft SUSY breaking parameters is to be welcomed, since it 
provides a window into the soft supersymmetry breaking Lagrangian. 
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APPENDIX A 
Notation and Conventions 



In this appendix we briefly summarize our conventions and explain in detail the 
notation concerning the diagonalization of the mass matrices. The superpotential of 
the MSSM+z/'^ model is given by : 

+yA{K)ABl%K - e\{\e)ABrBh^, + ^^Kh^a) + (54) 

where €12 = —£21 = 1, v4,i? = 1, ..,3 are flavour indices and a, (3 = 1,2 are SU{2)l 
indices. The Yukawa matrices are diagonalized by the following transformations : 

^"A„T"^ = a;, S'^XdT'^^ = X'a S^XeT^^ = K S^'X^T"^ = XI (55) 

where the primed A's are diagonal. In this notation the CKM matrix is given by 
V = 2""2"<^t Tj^g f^u lagrangian of the model also includes trilinear, soft scalar 
masses and gaugino masses given by : 

V = e^pi u'A{Au)ABqBhZ - d'A{Ad) ABQshd 

+D'},{A,)ABhhi - ~e^j,{A,)ABryd + m:,Kh^^) + h.c. (56) 

£ = l{ml\hu\^ + m\\hd\^ + q*A{ml)ABqB + l*A{m'{)ABh + 

UAimlc)ABUs + dA{rh%)ABdB + eAimlc)ABeB + ^a{^Ic)abVb) + h.c.(57) 

C = |AfiAiAi + iMsA^A^ + iMgA^A™ + h.c. (58) 

The chargino and neutralino mass matrices can be conveniently written in the basis 
of the following 4-component gaugino and higgsino fields : 
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B 



Thus, 



and, 



C 



„ , / M2 V2MwSp 



HI 



R 



HI 



+ h.c. 



(60) 



(61) 



/ Ml 



--(bw H,H^ 



-Mzcpse Mzspse \ ( B \ 

M2 Mzcpce -Mzspce 

MzCfsSg MzC/sCg -/i Hj 

-Mzspce -/i 



V Mzspse 



I 



(62) 



where sp = sin/3, C/3 = cos/3 (tan/3 = V2/V1 = the ratio of the up/down Higgs VEVs) 
and 6 the weak mixing angle. The chargino mass matrix M*^ in Eq. (^T]) and the 
neutrahno matrix in Eq. (|62D are diagonahzed by : 



Diag(m^-,m^- 



^ATj^Af^Aff _ Diag(m^o, m^o, m^o, m^o) 



(63) 
(64) 



The mass matrices for the charged scalar sparticles are written in the following basis : 

(65) 



dR 



Explicitly we find : 

(4 4 



m 



Ic + m„mi + M|Z„cC2/3 )\ur 



(66) 



(4 4) 

(el el) 



Tfi% + mdm\ + M|ZrfcC2/3 / V (ii? 
m^, + mem]; + M|ZecC2/3 / V e/? 



' mf + mime + MzZeC2f3 

-IJ,XeV2 + AeVi 

The light sneutrino mass matrix in the ul = u basis, after the heavy right-handed 
sneutrinos are integrated out, is given by : 



(67) 
(68) 



(mf + MfZ^^cs/j) Ul 



(69) 
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The Z factors in Eqs.(|66D-(|69D are defined by Zf = If — Qjsj where If is the isospin 
and Qf the electric charge of the / field : 

Zu = - {+l)sl Z^. = (0) - (-1).^ 
Z, = (-i) - {-l)sj Z,. = (0) - (+1).^ 

(70) 

Ze = (-1) - (-1).^ ^ec = (0) - ( + 1).^ 

Zu=m)-{ o)si 

The diagonalization of the up-type squark squared mass matrix M"^ of Eq. (|66|), 
down-type squark matrix M"'^ of Eq. (0), charged slepton matrix M^^ of Eq. ( |68D 



and sneutrino matrix of Eq. (|6^) is achieved in the following way : 

^-M^2^"t = Diag(mj^^,..,m^J (71) 

= Diag(msi,m£j,mt^,ms2,m52,mtJ (72) 

SV'S'^^ = Diag(m^^,..,m^J (73) 

= Diag(mj^,mg,,m^^,mj^,m,-2,m^J (74) 

S'M'^S'^ = Diag(m^^,..,m^J (75) 

= Diag(mei,m^^,m^^,me2,m^2,m^J (76) 

SIlMlISll = Diag(m^,,mi>2,mp3) (77) 



Finally it is convenient to define the following matrices that are, in a way, super- 
symmetric generalizations of the CKM matrix, i.e. they describe the fiavour properties 
of vertices that involve the interaction between a SUSY scalar particle and a standard 
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model fermion : 



{Vt)o.A 


= eL 


=1 '^aB+3'-'BA 


iVr')Aa 


= eL 


QU C"t 

=1 '-'AB'^B+Sa 


{vtU 


= eL 


cd rpd-\ 
-1 '-^aB-'-BA 


{Vi'')Aa 


= eL 


rpd cd^ 

=1 ^AB'-'Ba 


{Vt)aA 


= eL 


cd cd'\ 
=1 *-'aB+3^SA 




= eL 


cd cd] 
=1 ^AB*-'B+3a 


{Vt)aA 


= eL 


cu Td,\ 
=1 '-^aB-'-BA 




= eL 


T-id cut 
=1 -'-AB'-^Ba 


{VEDaA 


-EL 


--l{^LL)aBTBA 


{^LL)Aa 


= eL 


=1 ^IbI'S'l^bo 
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APPENDIX B 
The supersymmetry mass scale 



In this section we will define the effective SUSY scale Ms which we used to rep- 
resent the average sparticle masses. Above Ms the theory was described by the 
MSSM+i/^ and below it we had the Standard Model. Traditionally it is convenient 
to introduce three Tj scales that describe the effect of the decoupling of the SUSY 



particles on the gauge couplings p4 



rj, 11/25 9/25 4/25 1/25 /^^.x 

1 1 = rrij rrig (79) 



rr 9/25 8/25 4/25 3/25 1/25 /on\ 

J 2 = "^g f^T^ ^fj (°U) 



J 3 = m~ m 



1/2 1/2 



where rUg, mj, m^, ttlh, TUy^ and rrig are the squark, slepton, higgsino, heavy CP-even 
Higgs boson, wino and gluino masses respectively. In models with gauge unification 
at Mx the prediction for a<j(Mz) depends on Msusy given by : 

,^ Tn-25/19 TnlOO/19 Tn-56/19 /o^n 

Msusy = J 2 -^s \°^) 

Combining the equations above we obtain p| : 

3/19 / \ 3/19 / \ 4/19' 



/ \ 28/19 

Msusy = mJ'^] 



UK] ( 



(83) 



In this article we did not use the above scale to describe the average sparticle 
masses since as was emphasized in Refs. [|5|, Msusy is only related to the overall 
sparticle masses in the unlikely case of degenerate SUSY spectrum. Instead, we 
introduced a new scale Ms which was defined such that the sum of the squares of the 
threshold corrections induced by the decoupling of the SUSY spectrum on the gauge 
couplings is minimal : 



2 



^ ^E§lnf^)[ =0 where 5„ = (i,f,4) (84) 



dMs 2tt V Ti 
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Explicitly we find that Ms is given by : 

^ _ 2^225/1426 2^625/1426 2^576/1426 ^g^^ 
612/1426 228/1426 200/1426 156/1426 136/1426 34/1426 for\ 

= rriq rrig rrij (jib) 

where the exponent values of the ■mq,..,mH terms are : 0.43, 0.20, 0.14, 0.11, 0.10 and 
0.02 respectively. It is interesting to note that the unphysical nature of Msusy can be 
immediately identified through the appearance of negative exponents (see Eq. (|8^) 
whereas the Tj's in Eqs. ([79D-(pID and Ms in Eq. ( ^5]) and Eq. ( ^61) are weighed by 



positive numbers. 

In Table VI we list the low energy values of the masses of the gauginos, the 
effective SUSY scales in Eqs. ([791)- (pTD, and the value of the new supersymmetry 
scale Ms defined in Eq. ( |5B| ) against the often used Msusy of Eq. (P3|). We observe 
that while Ms is a good average of the Tj's, Msusy fails to represent a meaningful 
effective SUSY mass. 

TABLE VI. 



Case 






m-g 


Ti 


T2 


Ts 


Ms 


Msusy 


A 


173 


334 


942 


440 


454 


874 


588 


68 


B 


356 


674 


1787 


851 


887 


1650 


1132 


150 


C 


174 


334 


942 


438 


433 


879 


577 


53 


D 


356 


674 


1789 


818 


820 


1662 


1091 


103 


E 


173 


334 


976 


449 


474 


903 


610 


76 


F 


355 


673 


1845 


889 


935 


1700 


1181 


172 


G 


173 


334 


976 


451 


466 


907 


607 


68 


H 


356 


673 


1845 


876 


911 


1707 


1166 


150 


I 


264 


503 


1415 


718 


714 


1319 


915 


116 



TABLE VI. Predicted values for the bino (jrig), wino (jti^) and gluino (jrig) masses 
at low energies and for the three SUSY scales Ti and the effective Ms, Msusy scales 
defined in Eqs. (|86|),(^3|) (masses given in GeV units.). The first column indicates 
the input for each model as defined by list of values in Table I. 
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APPENDIX C 
Tables for the SUSY corrections 



In this appendix we present some examples which illustrate the magnitude of the 
SUSY corrections to the bottom and tau mass. 

In Table VII we systematically list the values of all the parameters, evaluated at 
low energy, that are needed to compute the gluino correction. The columns refer to : 
the input taken according to Table I; the value of the strong gauge coupling; the ratio 
of the MSSM VEVs; the mass of the gluino; the value of the Higgs mixing parameter 
fi] the masses of the physical bottom squarks; the value of the dimensionless function 
Ig = mgfil{m'^,m^ '"^ij''' ^ finally, in the last column, the value of the gluino 
correction (in percentage %.) In Table VIII we list similar values appropriate for the 
evaluation of the higgsino correction. For example, at = A^/47r, mi_^^ are stop masses 
and Ifj = Atfil{fi'^,m'^^,mf^). Comparing Table VII with Table VIII we read that 
A^mf, ~ —20% dominates over A^m^ ~ 6 %. 

In Table IX we present, in the last column, the values for the bino corrections to 
the tau mass A^rrir- In this table we also give the values for a' = g'^/An, for the bino 
mass m^, for the stau masses and for = m^fil{m'^^,mf_^,mlj. We read that, on 
average, A^rrir is -2.45 %. 

Finally we computed Ama and Arris which we show in Table X. Two comments 
deserve attention. Firstly, A^md^s — due to small Xd^s Yukawa couplings. Secondly, 
we see that the gluino contribution is not universal, i.e. A^m^ ~ A^uis 9^ A^mb 
due to mj- 9^ m^. In conjunction, they lead to non-universal SUSY corrections 
Anid^s 9^ Amf,, thus slightly affecting the ratio Xd,s{Mx) / Xb{Mx) ||69| . 

28 „, ^\ _ at/ \ix{x/y)+yz \n(y/z) + zx \ix{z/x) 

^ [X, y, Z) - {x~y){y^z){z~x) 
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TABLE VII. 



Case 


as 


tan/? 


m~g 


1^ 










A 


0.0904 


42.16 


942 


-417 


775 


691 


-0.305 


-24.65 


B 


0.0848 


44.16 


1787 


-795 


1431 


1316 


-0.311 


-24.72 


C 


0.0906 


35.94 


942 


-339 


790 


725 


-0.238 


-16.43 


D 


0.0850 


37.51 


1789 


-636 


1475 


1373 


-0.238 


-16.11 


E 


0.0960 


47.15 


976 


-494 


794 


692 


-0.358 


-34.36 


F 


0.0896 


49.40 


1845 


-942 


1452 


1323 


-0.367 


-34.44 


G 


0.0960 


41.99 


976 


-445 


804 


719 


-0.313 


-26.81 


H 


0.0897 


44.06 


1845 


-847 


1480 


1364 


-0.320 


-26.82 


I 


0.0920 


47.34 


1415 


-703 


1144 


1035 


-0.347 


-32.03 



TABLE VII. Values of the parameters required to compute the gluino correction to 
the bottom quark mass A^m^ ~ —20% (all evaluated at low energy.) The gluino 
mass is given by mg, the sbottom masses by mj^ ^ and the dimensionless function 
Ig is given by Ig = mgiJ.I{m1, m?^ , m?^). 



TABLE VIII. 



Case 


at 


tan P 


At 




mi. 


^2 






A 


0.0476 


42.16 


-780 


-417 


812 


665 


0.408 


6.51 


B 


0.0482 


44.16 


-1410 


-795 


1457 


1291 


0.401 


6.80 


C 


0.0363 


35.94 


-879 


-339 


831 


690 


0.385 


4.00 


D 


0.0364 


37.51 


-1605 


-636 


1504 


1345 


0.376 


4.09 


E 


0.0597 


47.15 


-728 


-494 


826 


676 


0.406 


9.10 


F 


0.0602 


49.40 


-1302 


-942 


1472 


1303 


0.399 


9.45 


G 


0.0499 


41.99 


-814 


-445 


842 


692 


0.416 


6.94 


H 


0.0504 


44.06 


-1465 


-847 


1506 


1338 


0.409 


7.23 


I 


0.0573 


47.34 


-1055 


-703 


1170 


1013 


0.400 


8.64 



TABLE VIII. Values of the parameters required to compute the higgsino correction 
to the bottom quark mass A^rub ^ 6 % (all evaluated at low energy.) The higgsino 
mass is approximately given by the stop masses are given by mj^ ^ and the 
dimensionless function 7^ is given by Ijj = At/i7(/i^,m|^,m|^). 
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TABLE IX. 



Case 


n' 


tan /? 


m f. 


/' 


m f ^ 


m f., 






A 


0.0104 


42.16 


173 


-417 


348 


178 


-0.669 


-2.33 


B 


0.0105 


44.16 


356 


-795 


653 


387 


-0.656 


-2.42 


C 


0.0104 


35.94 


174 


-339 


343 


207 


-0.507 


-1.51 


D 


0.0105 


37.51 


356 


-636 


657 


424 


-0.494 


-1.55 


E 


0.0104 


47.15 


173 


-494 


355 


143 


-0.860 


-3.36 


F 


0.0105 


49.40 


355 


-942 


652 


352 


-0.818 


-3.37 


G 


0.0104 


41.99 


173 


-445 


350 


175 


-0.715 


-2.49 


H 


0.0105 


44.06 


356 


-847 


655 


387 


-0.696 


-2.56 


I 


0.0104 


47.34 


264 


-703 


549 


329 


-0.641 


-2.52 



TABLE IX. Values of the parameters required to compute the bine correction to the 
tau lepton mass A^ttit- ~ —2.45%. (all evaluated at low energy.) The bino mass is 
given by m^, the stau masses by mf ^ j and the dimensionless function Ij^ is given 
by Ib = mstJ.I{m% , , m?^ ) 



TABLE X. 



Case 






A^nis 


A^TUs 


Arud 


Anis 


Am,}) 


A 


-19.24 


0.00 


-19.25 


0.01 


-19.24 


-19.23 


-17.20 


B 


-19.28 


0.00 


-19.28 


0.01 


-19.28 


-19.27 


-17.04 


C 


-12.95 


0.00 


-12.95 


0.01 


-12.95 


-12.94 


-11.35 


D 


-12.72 


0.00 


-12.72 


0.01 


-12.72 


-12.71 


-11.16 


E 


-26.38 


0.00 


-26.39 


0.02 


-26.38 


-26.37 


-24.40 


F 


-26.38 


0.00 


-26.38 


0.02 


-26.38 


-26.36 


-24.19 


G 


-20.92 


0.00 


-20.92 


0.01 


-20.92 


-20.91 


-18.87 


H 


-20.91 


0.00 


-20.92 


0.01 


-20.91 


-20.90 


-18.67 


I 


-24.40 


0.00 


-24.40 


0.02 


-24.40 


-24.38 


-22.55 



TABLE X. Gluino and higgsino SUSY corrections to the down and strange masses 
A^'^ma^s and the total correction to the down, strange and bottom quark masses 
Amdysfi (in percentage values.) 
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APPENDIX D 
Boundary conditions for the gauge and gaugino masses 



In this appendix we review the origin of Eq. (pO]) that relates the gauge and the 
gaugino masses of the SU{3)c ® SU{2)l ® f/(l)y MSSM and of the ^f/(4) ® SU{2)l ® 
SU{2)r Pati-Salam Model. 

We start by considering the constraint on the U{1)y, SU{2)f( and S'f/(4) gauge 
couplings g', g2R and g^. The covariant derivative for the heavy Higgs boson H that 
breaks the GUT symmetry is : 

D,H = d,H + ig^RT^RW^R^H + tg.T^G^H (87) 

where Tr = ^cr" and are the SU{2)ii and SU{A) group generators with associated 
gauge bosons and G"^ (a = 1, .., 3 and m = 1, .., 15). When H develops a non- 
vanishing VEV {H,^) = \/2V^ along the neutrino direction the quadratic interaction 
iDfj^Hl"^ generates the following mass terms : 

'^9lRK\r-Rr'R)22W^R''W'R^ (88) 

^92Rg4V,'{Tl,MT'%,W'R^Gf (89) 

2^2v;2(T'"T'^)44G""^G". (90) 



where r|, = diag(|, — |) and T^^ ~ y I |) |) ~|) are diagonal matrices. Upon 

explicit substitution of the group generators (and after adding to the above expres- 
sions similar terms associated with the H Higgs field, with a VEV {H^) = V2Vu) the 
mixing between the Wr, G^^ gauge bosons can be written as : 

where V'^ = + V^^. The matrix above can be diagonalized by a unitary matrix 
parameterised by the rotation angle a given by : 

92R vl5'4 . . 

sm a = I = cos a = = [92) 

lalR + lgl \l9lR + l9l 
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to yield a massless B state (the MSSM bino) and a heavy gauge boson X with a mass 
rrix = V'^{g2R + fs'l)- These are related to the W^, G^^ bosons through : 



B 
X 



cos a sin a \ / Wf^ 
sin a cos a ) \ 



G 



15 



COS a 
sin a 



sm a 



cos a J \X 



(93) 



Finally, in order to obtain the first equation in Eq. (|20|) we need to relate g' to (^2/? 
and (74. This can be achieved by examining, for example, the kinetic term for the 
fermions in the F representation : 



C = iFYD^F ~ -FY 94T^^Gf 



F ~ —wy^ 



1^4 (i) (sina5^)l M (94) 



and comparing it to the U{1)y neutral current Lagrangian : 



C = -uY 9' (I B, 



u 



We get : 



Thus combining Eq. 



and 



g = \l ^gi sm a 

we find : 

1 _ 1 1 

9'^ ~ qIr hi 



(95) 



(96) 



(97) 



which agrees with Eq. (|20|) after replacing g'"^ = ^gf. 

In the second part of this appendix we will justify the second equation in Eq. (pO|). 
We start by considering the following Higgs-gaugino-higgsino interaction : 



C = tV2H^ (^72Rr^A?e + g.T'^XT) H + h.c. 



(9^ 



where H is the heavy Higgs boson, H the heavy higgsino field, and A^, A™ the 
SU{2)ji, SU{4:) gauginos (a similar expression applies to H). When the Higgs bosons 
H, H develop their VEVs, we obtain (after appropriately substituting the generator 
matrices r|, T^^) : 



jO = -g2RH\%) WuH, + V,hJ\ + ^g^(-i\f) \v,H, + K^J + h.c 



(99) 
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Prom Eq. (|99D it is clear that the gauginos only couple to the VyHy+VyHy combination 
of higgsinos. Thus, it is convenient to define new fields A^^i and N2 through : 



N2 



cosOn sin6'7v 
— sin6'7v cos6'jv 



(100) 



-iAfj, —iXY', Ni,) basis can be written as : 



4 ? ^"Ij 





(101) 



where tan6'jv = V^/Vi,. The bosonic partners of A''i, N2 have VEVs given by {Ni) = 
V2V1 = y2Kcos^7v + v^Ksin^^ = V2{V^ + 1/2)1/2 ^^^^ ^jy^^ _ r^j^g gaugino 
mass matrix in the 



M4 

\92rVi -^/igiVi 
where the 13, 23 entries derive from Eq. (|99|) and M2R, M4 are explicit light soft 
masses for the SU{2)2r, S'f/(4) gauginos. The two heavy eigenvalues of the above 
matrix are approximately given by Mx^^ ~ ±{g2ji + ^glY^^Vi. The lightest eigenvalue 
(the bino mass) can easily be computed from Mi = DetMx/ (Mx^Mx^)- We find : 



g2RVi \ 
) 



Ml 



2^4 



gin + 1^4 ^ii? + 1^4 



rM. 



(102) 



Finally, replacing the result of Eq. (P7|) into Eq. ( |102[ ) gives us the expression we 
wanted to prove : 



g2R 'igi 

We would like to conclude this appendix with the following observation. The 
rotation that diagonalizes the W^-G^^ gauge bosons mass matrix (see Eq. (pi])) also 
simultaneously block diagonalizes the gaugino mass matrix in Eq. ( |101| ). Indeed we 
note that : 



M2R M4 

-.2 ~'~ 3^2 



/ 



V 



0' 

I M. 



A 



/ Mx, 0{M) 

0(M) 0(M) 0{Vi) 
V OiV^) 1 



(104) 



where = sin a, c„ = cos a and 0{M) is a number of order M2/j and/or M4. Thus 
we find from Eq. (|104|) that Mx^ = c'^M2r + s'^M^ which is nothing less than the bino 
mass of Eq. (|102|) re- written in terms of the gauge mixing angle a defined in Eq. (p^) . 
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APPENDIX E 
Expressions for the D-terms 



The D-term corrections to the soft masses in Eqs. (p6D -(^) arise when the rank 
of the Pati-Salam group is reduced from five to four due to gauge symmetry breaking. 
In this appendix we show that the coefficients of the corrections are related to the 
charge carried by the fields under the U{1) broken generator and that their magnitude 
depends on the difference between the soft masses of the heavy Higgs that break the 
GUT symmetry. 

We start by carefully reporting our index conventions. The matrix elements of 
the fields of the 5f/(4) ® SU{2)l O SU{2)r model in Eqs. (|15D-(|17D are indicated by : 
F^,: (4,1,2) (4,1,2) 

V: (1,2,2) (105) 

FT'^: (4,2,1) H'^'- : (4,1,2) 
where i,k,p,r = 1,..,4 are S'f/(4) indices, di,6,fj,a = 1,2 are SU{2)r indices and 
7,p = 1,2 are SU{2)l indices. The position of the indices is the following, the 
first /second index refers to the line/column of the matrices in Eqs. (p!5D-(p!7D. More- 
over up/down indices are related to the representation of the multiplet. For example, 
the lower p index of h indicates that h transforms in the 2 (anti-matter) representation 
of SU{2) under the SU{2)l symmetry, whereas the upper a index indicates that h 
transforms in the fundamental representation of SU{2) under the SU (2)^ symmetry. 

The D-term contributions from the SU{2)2r and S'f/(4) groups are given by : 

V = hln E D^,rD:^r^ + \gl E D^Dt. (106) 

a=l m=l 

We focus on the a = 3 and m = 15 contributions of Eq. ( |106D which involve the 



r| = diag(i, -i) T^' = diag(|, |, |, -i) (107) 

diagonal generators of the SU (2)2/?, and S'f/(4) groups. Using the notation of Eq. ( |105| ) 
we find that -D|/j, Dl^ are given by : 
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DIr = Hl{-4*)eAp + H^'^'ijD.^H^^ + F<^l(-rf + h]^ ^{rD^^h^ (108) 

Assuming that the heavy Higgs develop VEVs along the neutrino directions : 

{H) = = {H) = H^^ = (110) 

we can expand -Df^, D^' '■ 

DIr = Hli-rl;h,H. + Hli4)i,H, + 

d'\-r'R)ud^ + e^H-4niie' + u^\-rln,,u^ + v''\-rl*),^'^ + 

= - - l\df - ||gf + ||«f + i|z>f + l\h,\' - '^\hX (111) 

= ^AHW - - lldf + lief - ll^f + ||z>f + llgf - 1|/T](112) 

One can summarize the results of Eqs. (|111|) -( |TT^ ) by writing : 

dIr = Dh + y: im' Dl' = ^,Dh + /IE (113) 

where Dh = i(|i7,^p — ji^iyp) and denotes any of the light fields u'^, d'^, e'^, i)'^, q, 
I, hu, hd- The factor refers to the charge carried by with respect to the SU {2)r 
group and {B — L)/2 to the semi- difference between the baryon and lepton numbers 
of (f). These are can be read from the coefficients of the terms in Eqs. (|lllD -( |ri^ ) and 
are collected in Table XI. 
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TABLE XL 



B-L 



u'^ d!^ v'^ q I hu hd 

+ 1 _i -1 +i -i +i 

^2 2 2^2 2^2 

_i _i +1 +1 +1 -i 



TABLE XL Charges carried by the hght scalar fields under the SU {2)2r and (Baryon- 
Lepton)/2 symmetries. 



Finally, using Eq. ( p.l3| ) into Eq. (|106|) gives : 



V = IglADln + ^glDl'Dl' (114) 

= '^glRl'^DHT.um + bll'^lDHT.i^)^] (115) 

= km'~W)T.{9lRl^ + l9l{^)A\<P\' (116) 

The equation above deserves two comments. Firstly, we see that the broken U{1) 
generator X resulting from the GUT symmetry breaking : 

SU{4) ® SU{2)2B. ^ SU{3)c ® U{1)y ® U{l)x (117) 

is given hj X = I + {B — L) /2 whereas the unbroken hypercharge is F = — / + {B — 
L)/2. Secondly, comparing Eq. (|116|) with Eqs. (|26|)-(|33D we find : 



D' = 1{\H,\' - (118) 

From Eq. ( |118|) alone, one might be tempted to conclude that the natural scale for 
the D-term is of the order of the mass of the heavy Higgs. However this is not true. 

The scale of the D-term in Eq. ( |ll(j| ) can be estimated upon the minimization of 
the heavy Higgs potential given by : 

V = U9m + -29l){Hl - Hlf + \1{H,H, - Mlf + m\Hl + mlHl (119) 
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where the first term is a D-term, the second an F^-term where S is the gauge singlet 
of Eq. ( jlSD and Mh ~ 10^^ GeV, and m^j, mn soft scalar masses. 

The minimization conditions for V, re- written as combinations of dV/dS^ ± 
dV jdHy = 0, are given by : 

[gl{Hl - Hlf + \1{H,H, - Ml)]{H, + H,) = -mlH, - mlH, (120) 

[gUHl + Hlf~\l{H,H,-Ml)]{H,-H,) = -m\H, +mlH, (121) 

where gjj = \{g2R + ^gf)- In the hmit of neghgible soft masses the right-hand side of 
Eqs. ([T20|)-([T2lD vanishes. Thus, we find that = = Mh, implying that D = 0. 
This result could had been anticipated since, if mjj = rriH then the Higgs potential 
is invariant under the exchange of the ^ H,j fields. When the soft masses are not 
zero, but still much smaller than Mh a perturbative solution to Eqs. (|120|) -( p!2l| ) is 
appropriate. In terms of the new parameters m, fh defined by : 

H^ = Mn-rh H^ = Mh -m (122) 

the minimization conditions are : 

2Xl{Th + m)M = mff + m% Ag]j{m - m)M = m\ - m]j (123) 

Thus we find : 

1 ( mW 1 f mlrX 



and finally : 



2 mjj — rrig 

^ = 4A| + 2gl,, + 3gl ^^^^^ 



We see that, in spite of D being the difference of two GUT scale masses, it actually 
scales with the difference between the heavy soft Higgs masses. The reason is because 
~ eM|f where e = {mj^ - m^) /Mfj is a very small parameter which measures the 
amount of the H,y Hy symmetry breaking of the potential V in Eq. ( |119|) . 
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